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ABSTRACT 

Complementary one-, two-, and three-dimensional tests for detecting the presence of 
substructure in clusters of galaxies are applied to recently obtained data from the 
2dF Galaxy Redshift Survey. The sample of 25 clusters used in this study includes 16 
clusters not previously investigated for substructure. Substructure is detected at or 
greater than the 99% CL level in at least one test for 21 of the 25 clusters studied here. 
From the results, it appears that low richness clusters commonly contain subclusters 
participating in mergers. About half of the clusters have two or more components 
within 0.5 hr l Mpc of the cluster centroid, and at least three clusters (Abell 1139, Abell 
1663, and Abell S333) exhibit velocity-position characteristics consistent with the 
presence of possible cluster rotation, shear, or infall dynamics. The geometry of certain 
features is consistent with influence by the host supercluster environments. In general, 
our results support the hypothesis that low richness clusters relax to structureless 
equilibrium states on very long dynamical time scales (if at all). 

Key words: galaxies: clusters: general - large-scale structure of Universe 
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1 INTRODUCTION 

It is well recognized that many clusters of galaxies ex- 
hibit substructure where the operational definition of 
substructure adopted here will be departures from oth- 
erwise smooth galaxy density and velocity equilibrium 
distributions. Representative previous st udies of samples 
of several clusters are given by, e . g.. IWest & Bothun 
Jl99d).lRhee. van Haarlem, fc Katgertl Jl99ll) . iBirdl | l994" 
lEscalera et"all (Il994l) . IWest. Jones fc Formar] Il99. 
Pinkney, Roett iger, Burns, and Bird (1996), |Gjrajdj^et_alJ 
(1997), and ISolanes. Salvador-Sole, fc Gonzalez-Casadc 
(1999). 1 The presence of subclusters affects estimates of 



quantities such as the dynamical mass and mass-to-light 
ratio, the mean gravitational potential, and the global 
ellipticity. Substructure may also occur in the intracluster 
gas distribution traced by X-ray emission, and cross- 
correlating galaxy-gas structure is an i mportant part of 
characterizing cluster dyna mics (see, e.g.. lDavis et aljll995l 
iBird. Davis, fc Beers! Il995l or references contained in the 
review article bv iRosati. Borgani. fc Normanl 120021) . In 
addition to determining the dynamical state of a cluster 
at a single epoch, the interaction of subclusters within the 
cluster (including mergers) can significantly affect dynami- 
cal evolution. Identifying and quantifying substructure may 
also reveal clues as to the initial conditions of the density 
perturbations from which clusters form and evolve. 

The purpose of the present study is to apply a com- 
prehensive, yet reasonably small, set of tests for substruc- 
ture to recently obtai ned data from the 2dF Galaxy Red- 
shift Survey (2dFGR S) JCollesJl99SUColless et all200ll and 
iMaddox et al.l [*1998). With the large number of redshifts 
provided by the Survey, it is possible to explicitly correlate 
the substructure detected in the projected two-dimensional 
surface distribution with that found from three-dimensional 
tests combining positional and velocity information. This is 
usually sufficient to distinguish true substructure from pro- 
jection effects including foreground and background contam- 
ination. Thus, the work presented in this paper is intended 
to provide survey results, identify potential trends and com- 
mon characteristics, and highlight features warranting fur- 
ther study; it is not intended to necessarily provide a defini- 
tive last word. 

Among its goals, the 2dFGRS is designed to allow in- 
vestigation of the properties of galaxy groups and clus- 
ters by providing a large, homogeneous sample in redshift 
space. The source catalogue used as the Survey base is es- 
sentially a revised and extended version of the APM cat- 
alogue with target galaxies of extinction-corrected magni- 
tudes bj < 19.45 and redshifts with a m easured RMS un- 
certainty of 85km s _1 iColless et al .1200 J) . The median red- 
shift of the catalogue is z — 0.11. A preliminary study of 
the Survey database identified 431 Abell, 173 APM, and 
343 EDCC clusters and provides p recise redshifts, velocit y 
dispersions, and updated centroids iDe Propris et alJ [2002). 
Additionally, this new data affords an opportunity to as- 
sess the completeness of, and contamination in, the older 



1 Because of the large number of published substructure studies, 
it is not possible to cite all of them; the references herein are 
meant to be representative and directly relevant to the results 
presented in this paper. 



catalogues. The 25 clusters of galaxies listed in Table to 
be studied here are selected from this new compilation, and 
include 16 clusters not previously investigated for substruc- 
ture. The proper distances to each cluster at the epochs of 
emission and observation, d p (t e ) and d p (to), are calculated 
under the assumption of a flat universe with S1a =0.7 and 
fl m = 0.3. All distances and linear scales are given in units 
of h' 1 Mpc where h = H /W0. 

The last two columns of Table contain the original 
richness classification (where available) as well as a revised 
richness based on the current 2dFGRS catalogue of clus- 
ters. Based on the memberships in the present catalogue, 
it can be seen that under the st rict definitio n of the rich- 
ness classes proposed by Abell iAbelll Il95fil) . 9 of the 25 
clusters in this study would be classified as "large groups" 
instead of low richness clusters. Nevertheless, these struc- 
tures are included in this study under the (loose) heading 
of "low richness clusters". Also note that since the mem- 
bership completeness here is estimated at 80 — 90%, it is 
possible that 100% completeness would alter these richness 
estimates. 

Among the previous substructure investigations in- 
vo lving relatively large samples of clusters is the study 
bv lSolanes. Salvador-Sole, fc Gonzalez-Casadol lll999T) . This 
work considered subclustering in 67 rich clusters contained 
in the ESO Nearby Abell Cluster Survey (ENACS) cata- 
logue. Using the best data available at the time, the ENACS 
catalogue is a reasonably homogeneous dataset with a com- 
pleteness of ~ 60 — 80% (fraction of cluster members with 
redshifts), less but similar to that of the 25 2dFGRS clus- 
ters selected for this study. The completeness for the sample 
here is estimated to be ~ 80 — 90% for each cluster, and the 
sample contains all clusters in the 2dFGRS catalogue with 
N > 70. Two important differences between the ENACS 
study and this work are (1) the different test suites used to 
test for substructure, and (2) the average number of cluster 
members in the ENACS study is significantly less than here: 
the average membership in the ENACS study was N = 45 
with 48 of the 67 clusters having N < 50, whereas the av- 
erage membership number for this study is N — 99 with 
the minimum number being N = 73. This is important due 
to the dependence of the reliability and sensitivity of the 
statistical tests as a function of N. 

The paper is organized as follows. Section 2 presents 
the method by which the ellipticity parameters, core radius, 
and central density are estimated. The tests used to detect 
the presence of substructure are described in Section 3 with 
a detailed analysis of the individual clusters presented in 
Section 4. A summary of the results and a concluding dis- 
cussion are given in Section 5. Finally, Appendix 1 contains 
visualization plots for each cluster. 



2 GLOBAL CLUSTER PROPERTIES 

Prior to applying tests for substructure, it is necessary to 
characterize each cluster in terms of global parameters such 
as a core radius and maximum central density, mean velocity 
and velocity dispersion, and global ellipticity and position 
angle. Here, "global" refers to the calculation of properties 
over an angular area that contains most or all cluster mem- 
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Table 1. The sample of 2dFGRS clusters selected for substructure analysis. Where avail- 
able, and using the original defini tions given by Abell, the origin al ACO richness classifi- 
cation R, is given iAbelll958l and lAbell. Corwin, fc Olowirll989t) along with the richness 
class determined from the 2dFGRS catalogue of clusters. The number in parentheses 
in the R2dF column is the number of galaxies between m.3 and + 2 in the current 
catalogue. The proper distances d p (t e ) and d p (to) are given in h~ 1 Mpc. 
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bers. First, the cluster selection and foreground/background 
rejection methods are briefly summarized. 



Essentially, this first catalogue of 2dFGRS clusters is a 
preliminary catalogue based on comparing the 2dFGRS re- 
sults with previously identified c lusters of ga laxies sourced 
from the catalogues of Abell j ^bejlj and the re - 

vised /s upplemented Abell llAbell. Corwin. fc 01o winlll989l). 
APM dPalton et alJ Il997l) . and EDCC jLumsden et all 
1992). The 2dFGRS catalogue was searched to identify clus- 
ters with centroids given in the above catalogues within 1° of 
the centre of an observed survey tile. If the centroid of a cat- 
alogued cluster was found in a 2dFGRS tile, the 2dFGRS 
redshift catalogue was then searched for objects within a 
specified radius of the centroid. This process isolates a cone 
in redshift space containing candidate cluster members as 
well as foreground and background galaxies. The list of can- 
didate members was then refined by inspection of the Palo- 
mar Observatory Sky Survey (POSS) plates. Finally, fore- 
ground and background contamination was eliminated (or 
significantly reduced) through subsequent detailed analysis 
of the redshift cone diagrams and, when necessary, redshift 
histograms. Further details concern ing cluster selection can 
be found in lDe Propris et alJ (I2002I) . 



2.1 Cluster Centroid and Ellipticity 

Each galaxy position in (RA, Dec) is converted to cartesian 
(x, y) coordinates with relative separations calculated using 
the proper distance at photon time of emission, d p (t e ). The 
cluster centroid is determined using an iterative search algo- 
rithm. A trial centre is computed from the arithmetic mean 
of all galaxies in the cluster, and is then recalculated using 
only the galaxies within an approximately 0.5 h~ Mpc ra- 
dius of the initial centre. This process is repeated until it 
converges and the final centre determined. The reason for 
restricting the centroid computation in this manner is to 
avoid bias from large substructures well removed from the 
central regions. 

For a quantitative measurement of cluster shape, it 
is convenient to use the dispersion ellipse of the bivari- 
ate normal frequency fu nction of position vectors (see 
iTrumpler fc We aver! Il9 52|) and f irst u sed to study cluster 
ellipticity bv ICarter fc Metcalfe! (|l980). The dispersion el- 
lipse is defined as the contour at which the density is 0.61 
times the maximum density of a set of points distributed 
normally with respect to two correlated variables although 
the assumption of normally distributed coordinates is not 
necessary to achieve accurate shape parameters. Using the 
moments 

M10 = x (la) 
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Table 2. Cluster Ellipticity and Position Angles. Position angles are measured counter- 
clockwise from north (positive y-axis). 
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To achieve the best results, an iterative approach is 
adopted here similar to th at used in previou s studies such as 
ICarter fc Metcalfe! (ll98CT> or lBurgetti lll982l) . An initial circle 
of a given radius R is defined about a trial centre and a new 
centre and dispersion ellipse are then calculated using only 
the galaxies contained within the circle. The semi-major and 
semi-minor axes A and B of the new ellipse are constrained 
to satisfy R = V AB. The process is repeated until con- 
vergence is achieved (usually requiring less than four iter- 
ations). It is possible for the iterative solution to become 
trapped in a local extremum, but this can be avoided by 
tracking each iteration. In applying the algorithm to the 18 
richest clusters in the Dressier 1980 catalogue, it was noted 
that the result s are sensitive to the presence of substructure 
jBurgetdll982T) . The ellipticities and position angles for the 
25 clusters are shown in Table [5] for four different average 
distances from the cluster centroids; the dispersion ellipses 
for mean radii of 0.5 and 1.0 h _1 Mpc are superimposed on 
the plots of galaxy positions in Appendix A for each cluster 
analyzed in Section 4. 



Since a model density distribution can be significantly 
affected by ellipticity, this must be accounted for when con- 
structing Monte Carlo catalogues of random clusters for 
comparison against actual data. Also note that it is some- 
times nontrivial to distinguish true ellipticity from subclus- 
ter bias, and, in fact, the algorithm used to compute ellip- 
ticity can be used to also probe for substructure. 
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The position angle of the major axis with respect to north 
is then 



Mil / 2 



(3) 



with Fa (> Fb) being the semi- major axis, and the elliptic- 
ity is 

F B 



e= 1 - 



r i 



(4) 
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Table 3. Cluster Core Radii and Central Densities. The errors quoted are for the par- 
ticular fit. 
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2.2 Density Profiles and Core Fitting 

The assumed models for the density and velocity distri- 
butions are those for a cluster in dynamical and thermal 
equilibrium, viz., the empirical King approximation to an 
isothermal sphere and a Gaussian velocity distribution. It 
should be remembered that while the projected number den- 
sity distributions of many clusters are fit well by King pro- 
files, the core radius and central density are fitting variables 
that are convenient, but arbitrary, parameters for character- 
ising the compactness of cluster central regions. In particu- 
lar, the core radius has no intrinsic dynamical significance. 
Thus, the results presented in Table|H]should be regarded as 
comparative geometrical measures across the sample with- 
out dynamical implications. Similar considerations apply to 
the usage of the term 'core' throughout this paper. 

For a circularly symmetric c luster proje ction, the two- 
dimensional surface King profile dKindll9'721) takes the form 



where A and B are the semi-major and semi-minor axes, 
respectively, and where it is conventional to define an equiv- 
alent mean core radius, R c — V AB. 

As a starting point, values of R c and ao for the circu- 
lar profile of Eqn. J5J are calculated using a least squares 
method. First, the expected galaxy count n(r) at a given 
distance from the centre is found by integrating Eqn. |J5J, 

n(r) = irR 2 c a ln[(R 2 c +r 2 )/R 2 c ] , (7) 

and inverting for the radius r(n) of the circle containing n 
galaxies, 



r(n) = R c \/exp (n/ivR^ao) — 1 . (8) 
Each galaxy then becomes a data point in the merit function 



x 2 (R c ,a )=J2[- 



N , ... N 2 

n — r(i) 



(9) 



<r(r) 



(To 



1 + 



(5) 



where a(r) is the galaxy density at a distance r from the cen- 
tre of the cluster, ao is the maximum galaxy density of the 
cluster, and R c is the core radius of cluster. By extension, 
the two-dimensional King profile for an elliptically symmet- 
ric surface density distribution is 



a{x,y) 



<T() 



t /+(sr 



i+ 4 



(6) 



where is the actual distance of the ith galaxy from the 
centre of the cluster and r(i) is given by Eqn. The min- 
imisation of the merit function requires the derivatives with 
respect to parameters R c and ao to be set equal to zero, 



dR c 



(0 



dr 



dR c 







and 
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where 
dr 
~dR c 

and 



(ivaoRc — n) exp (n/ivR^ao) — ivR^ao 
ira R 2 \/exp (n/TvR 2 a ) — 1 



dr —nexp(n/-KR^ao) — ivR^ao 
da 2na R? c v'exp (n/irR 2 a ) - 1 



(11a) 



(lib) 



The Levenberg-Marquardt algorith m is t hen used to find the 
best fit for R c and cr l|Press et al.lll992h . 

The circular fit above requires modification when the 
cluster core is highly elongated or multimodal. While it is 
straightforward to extend the fitting method above for the 
four free parameters in the elliptical King profile of Eqn. ©, 
it is possible to employ a simpler approach. If an ellipse with 
ellipticity e is constructed such that V 'AB = lh^ 1 Mpc = 
V Ar-iBr-i = R=i, then the coordinates of the perimeter 
satisfy 



A 2 



+ 



B 2 

D R=1 



(12) 



Given the coordinates (xi, yi) of the ith galaxy in the cluster, 
there exists a concentric ellipse which includes the coordi- 
nate (xi,yi) defined by 

1 (13) 



a) + Et 



This ellipse is defined by = e and Ri — V ' AiBi with Ri 
satisfying 



ft = R=i 



A 2 

A R=1 



+ 



R 2 



(14) 



The circular core fitting technique described above can 
now be utilized in the following manner: 

(i) Calculate values for the semi-major axis Ar = \ and 
the semi-minor axis -Bh=i using values for the ellipticity e 
from Table |5] (this requires a certain amount of judgement 
to select the best values), 

(ii) Transform the (xi,yi) coordinates for each galaxy to 
(x'i,yi) for the coordinate system rotated to have the :r'-axis 
along the semi-major axis of the ellipse, 

(iii) Calculate for all galaxies the quantity 



Ri — R- 



=i + M) 2 /b 2 r __ 



(15) 



(iv) Sort the values of Ri in ascending order. This new 
order of galaxies, denoted by n, along with their values for 
Ri(n), gives the number of galaxies found on or within the 
ellipse where R — Ri(n), 

(v) Use the above values for Ri(n) in place of r(n) used 
by the circular fitting algorithm to obtain new values for R c 
and cro from the least-squares fit, 

(vi) Recompute A c and B c using the original e = 1 — B/A 
and R c — V ' A c B c . 

The results for elliptical surface density profile fits are shown 
in Table|3]for the 25 clusters. Comparisons of the results for 
the circular and elliptical fits shows that the fit to a circular 
profile provides a fair representation of the number distri- 
bution even when the actual shape is noticeably elongated. 
On average, compared to the elliptical profile, the circular 
fit also is better over the extent of the entire cluster. This 



can be due to the cluster ellipticity generally decreasing as 
a function of distance from the centroid as well as the circu- 
lar fit essentially averaging over clumpy distributions that 
might bias the elliptical fit. However, the elliptical fit repro- 
duces the density profile in the central regions significantly 
better than the circular fit for clusters such as, e.g., Abell 
4012 (Figure Q). 

For some clusters containing multi-component core re- 
gions such as Abell 930, the circular fit is slightly better than 
the elliptical fit in the central regions. This is presumably 
due to both components containing approximately the same 
number of galaxies so that the computed centroid is not at 
the centre of either component. Further detailed co mpar- 
isons and fitting plots for all clusters can be found in IVickl 
(2004). 



3 TESTS FOR CLUSTER SUBSTRUCTURE 

The dangers of inferring the presence of substructure based 
on visual inspection of contour or gray scale plots alone is 
well known, yet this remains an appropriate component of a 
comprehensive test suite in order to better guide and inter- 
pret a statistical analysis. Two early studies utilizing contour 
plots to detect possibl e substructure i n clus ters of galax- 
ies we re conducted by iGeller fc Beers! (Il982ft and iBurget tl 
<1982f). In c r eating contour plots of clusters drawn from 
the iDresslerl lll980l) catalogue, the Geller and Beers study 
applied a "boxcar smoothing" technique that was a best- 
available technique at that time. The improvements gained 
using an adaptive kernel approach to co ntour plotting are 
shown for many of these same clusters in lKriessler fc Beersl 
il997f) . In addition to number density contours, the Bur- 
gett study presented luminosity-weighted contour plots of 
clusters in an attempt to correlate light with mass. As sub- 
structure studies have evolved with experience and numer- 
ical sophistication, it is now conventionally accepted that 
a synergistic combination of quantitative tests is necessary 
to detect different types of structure as well as to mitigate 
the results of false positives occasionally given by almost all 
estimators. 

The test suite adopted for this study consists of a 
variety of one-, two-, and three-dimensional tests. These 
include visualization plots such as contour and nearest 
neigh bor plots, velocity d istributio n statistics, and the a 
test JWest fc Bothunlll990l). 8 test JWest. Oemler. fc Dekel 
Il988h . and k test iColless fc Dunnll996t) . For the latter three 
statistical tests, and for all clusters, a sample of 10,000 sim- 
ulations was used to normalize the test and to determine the 
statistical significance. Use of the two point angular correla- 
tion function and the Fourier elongation test were considered 
but rejected on the grounds that they would not provide any 
new information compared to the chosen test suite. 



3.1 Velocity Statistics as a Substructure Indicator 

The velocities for the galaxies in a relaxed cluster in thermal 
equilibrium are expected to be distributed normally with re- 
spect to any one cartesian coordinate. However, it should be 
recognized that taken by itself, an apparently non-Gaussian 
velocity distribution may reflect either a cluster not in equi- 
librium or a cluster containing one or more dynamically- 



Substructure in 2dFGRS Clusters 




Figure 1. Comparison of the circular and elliptical King fits for Abcll 4012. The solid 
line represents the expected count n of galaxies for the fitted density model within cither 
a circle of radius r(n) (circular King profile) or a circle of geometric mean radius R(n) = 
[A(n) £?(«)] 1//2 for an ellipse with semimajor and semiminor axes A(n), B(n) (elliptical 
King profile). The points in each plot represent a galaxy at a specified distance from the 
cluster ccntroid. Points above the fit line are outside the test circle, but should have been 
inside at that distance; points underneath the fit line are within the test circle, but yield 
a higher density than the model density at that distance. 



Table 4. Statistical moments of the cluster velocity distributions, equivalent standard 
deviations for the skewnness and kurtosis, and the confidence level (CL) derived from 
the Xv e i result for consistency with a Gaussian distribution. The entry for Abell S333* 
presents the velocity statistics for Abell S333 without the 8 galaxies in the group having 
cz < 18,500 km/s, and the entry for Abell 4013* presents the velocity statistics without 
the 14 galaxies having cz > 17, 500 km/s. 



Cluster 




Skewness 


#of 


Kurtosis 


#of 


CL to reject 






(km/s) 




o\s 






Gaussian (%) 


Abell 


930 


856 


0.30 


0.7 


-1.00 


1.0 


>99.9 


Abell 


957 


708 


0.31 


0.8 


-0.45 


0.4 


49 


Abcll 


1139 


484 


-0.21 


0.6 


0.18 


0.2 


76 


Abell 


1238 


551 


0.30 


0.7 


0.04 


< 0.1 


14 


Abcll 


1620 


1001 


0.77 


1.9 


0.65 


0.6 


97 


Abcll 


1663 


881 


-0.60 


1.5 


0.86 


0.8 


36 


Abcll 


1750 


897 


0.43 


1.0 


-0.87 


0.8 


98 


Abell 


2734 


984 


0.71 


2.1 


1.38 


1.6 


85 


Abell 


2814 


894 


-0.14 


0.3 


0.09 


0.1 


68 


Abell 


3027 


838 


-0.57 


1.4 


-0.76 


0.7 


99.9 


Abell 


3094 


728 


0.42 


1.1 


-0.33 


0.4 


93 


Abcll 


3880 


784 


-0.33 


0.9 


-0.02 


< 0.1 


36 


Abcll 


4012 


471 


-0.08 


0.2 


0.37 


0.3 


96 


Abell 


4013 


854 


0.93 


2.2 


-0.12 


0.1 


>99.9 


Abcll 


4013* 


455 


0.07 


0.2 


0.15 


0.1 


97 


Abell 


4038 


835 


-0.04 


0.1 


-0.31 


0.4 


96 


Abell 


S141 


403 


0.82 


2.2 


0.58 


0.6 


98 


Abell 


S258 


557 


0.07 


0.2 


-0.52 


0.5 


70 


Abell 


S301 


679 


-0.01 


< 0.1 


1.92 


1.9 


99.9 


Abcll 


S333 


933 


-1.13 


2.5 


1.04 


0.9 


99.9 


Abell 


S333* 


573 


0.01 


< 0.1 


-0.93 


0.8 


93 


Abell 


S1043 


1271 


0.86 


2.3 


-0.53 


0.6 


>99.9 


APM 


268 


769 


-0.65 


1.7 


-0.47 


0.5 


98 


APM 


917 


471 


0.18 


0.4 


0.03 


< 0.1 


73 


APM 


933 


1033 


0.77 


2.2 


0.81 


0.9 


97 


EDCC 


365 


531 


0.72 


1.7 


1.00 


0.9 


83 


EDCC 


442 


726 


0.00 


< 0.1 


1.72 


2.1 


51 
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Figure 2. The skewness values for the 25 clusters with the stan- 
dard deviation of the computed skewness value for a single cluster 
being approximately y/l5/N. From the figures it is evident that 
there is an excess of high skewness values over that expected if 
they were normally distributed. 



6 



24 
20 
16 



JD 

E 12 



-10 12 

Kurtosis - 3 



1 2 3 4 5 
I Kurtosis-3 I /V(96/N) 



Figure 3. The kurtosis values for the 25 clusters with the stan- 
dard deviation of the computed kurtosis value for a single cluster 
being approximately From the figures it is evident that 

the velocity distributions for the majority cannot be distinguished 
from Gaussian distributions based solely on kurtosis values. 



bound subgroups embedded in an otherwise smooth, approx- 
imately equilibrium, distribution. Since all semi- invariants 
of higher order than two vanish identically for a Gaussian 
(normal) distribution, deviations from normality can be es- 
timated by calculating the skewness and kurtosis (quantities 
proportional to the third and fourth order semi-invariants). 
In the following, all cluster redshifts have been transformed 
to velocities via the standard relativistic Doppler shift for- 
mula. 

The skewness is given by 



1 



(16) 



with v and a the mean velocity and standard deviation de- 
termined from the observed line-of-sight velocities Vi of the 
iV cluster members. A positive (negative) value of S implies 
the distribution is skewed toward values greater (less) than 
the mean. The kurtosis coefficient is defined as 



K 



i N 

w5> 



Because the kurtosis of a normal distribution is identically 
equal to 3, kurtosis values are conventionally presented 
with this subtracted from the result. Note that positive 
(negative) values of K indicate distributions more strongly 
peaked (flatter) compared to normal distributions. Figures 




[5] and [3| show the distribution of the skewness and kurtosis 
values for the 25 clusters as well as estimates of their 
statistical significance in terms of the number of standard 
deviations away from the Gaussian values. Examples of 
previous studies of cluster structure using computed skew- 
ness a nd kurtosis to characterize the velocity distr i butio n 
mcludelBird fc Beers! dl993h. I Ashman Bird, fc Zepj <1994h 
' "1994). I West. Jones fc Format] dl995l). iPinknev et al] 
and Solancs, Salv ador-Sole, fc Gonzalez-Casada 

In practice, for a cluster containing N ~ 50— 100 mem- 
bers, the skewness and kurtosis can provide suggestive but 
not necessarily definitive estimates as to whether a given dis- 
tribution is significantly non-Gaussian. In order to obtain a 
better quantitative check for normality, it is convenient to 
use a \ 2 statistic to estimate the probability that a test dis- 
tribution can be distinguished from a Gaussian with the null 
hypothesis that the test distribution is, in fact, normal. To 
be reliable, the x 2 algorithm requires the binned test distri- 
bution to have a minimum number of datapoints in each bin 
to guarantee a reasonable probability for occupancy. This 
can be accomplished either by shifting/changing the num- 
ber of fixed-width bins, or by using variable- width bins with 
(chosen) occupation probabilities above the minimum ac- 
ceptable level. In testing the velocity distributions here, the 
adaptive method was implemented with a bin width varied 
to yield a minimum 5% occupation probability for each bin. 

Table [I] shows the velocity dispersion a ve i, skewness 5", 
the number of standard deviations as away from the skew- 
ness value expected for a normal distribution, kurtosis K, 
the number of standard deviations <jk away from the kurto- 
sis value expected for a normal distribution, and the con- 
fidence level CL derived from the Xvei test to reject the 
hypothesis of a normal distribution. Inspection of the re- 
sults shows that 13 of the 25 clusters have catalogue distri- 
butions inconsistent with Gaussian distributions at > 95% 
confidence level (CL). Comparison of the Xvei statistic with 
the skewness and kurtosis values shows no consistent cor- 
relation as to whether a distribution is (or is not) normal. 
In general, the skewness and kurtosis are not necessarily 
reliable indicators of departures from normality for these 
cluster membership sizes. However, it should be noted that 
even the x 2 test is not 100% reliable, and the results depend 
somewhat on the degrees of freedom used in each case. 



3.2 Contour, Segmentation, and Nearest 
Neighbor Visualization Plots 

In spite of the dangers of giving too much weight to visual 
appearances, contour and other visualization plots remain 
useful when attempting to correlate the results of purely sta- 
tistical tests to spatial structure or when comparing galaxy 
structure to gas structure revealed through X-ray emission. 
Thus, in the analysis of the clusters considered in this study, 
four different types of two-dimensional visualization plots 
are included in addition to the scatter plot of the galaxy 
positions. 

The number density contour plots for the clusters 
are constructed with galaxy positions binned in 0.3 x 
0.3 (/i _1 Mpc) 2 square cells and then smoothing the con- 
tours with a bicubic interpolation algorithm. In order to 
provide a sufficient number of contour lines to give a visual 
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sense of cluster structure, it was found that the sample could 
be broken into two general categories depending on whether, 
for a given cluster, any 0.3 x 0.3 (hT 1 Mpc) 2 cell exceeded a 
count of 8 or more galaxies. For these clusters, contour levels 
were drawn at 20 gal/(ft~ 1 Mpc) 2 intervals, while for clus- 
ters with no cell exceeding this threshold, contour lines were 
drawn at 10 gal/(/i _1 Mpc) 2 intervals. In the first category 
(20 gal/(h -1 Mpc) 2 ) are Abell clusters 930, 957, 1139, 2734, 
1663, 2734, 2814, 3094, 3880, 4012, 4013, 4038, S141, S258, 
S301, S1043, APM clusters 268, 917 and 933, and EDCC 
442. In the second category (10 gal/(/i _1 Mpc) 2 ) are Abell 
clusters 1238, 1620, 1750, 3027, and S333 and EDCC365. 

As a means of comparing the cluster luminosity dis- 
tribution to the number density distribution, we have also 
constructed luminosity-weighted contour plots. As for the 
number density contour plots, the initial cell size is 0.3 x 
0.3 {hT 1 Mpc) 2 . The fc-corrected absolute magnitude for 
each galaxy, Mj, is calculated using th e standard T aylor 
expansion to first order in z (see, e.g., |Peebleslll993l pp. 
328-330), 

Mj = bj-k- 42.386 + 5 log 10 h - 5 log 10 z 

-1.086(1- go)* + ... , (18) 

where qo — —0.55 is the deceleration parameter value for 
f2 m = 0.3 and S1a = 0.7. The fc-corr ections are com- 
puted from the emprical fits presented in lDe Propris et alJ 
(2004). Each contour line represents an increase of ~ 5 x 
10 10 L o /(/i~ 1 Mpc) 2 (except for A2814 and A1750 which 
were set to ~ 1 x 10 11 LqUJtC 1 Mpc) 2 ), and, as for the 
number density contours, the luminosity density contours 
are smoothed using bicubic interpolation. 

The segmentation plots included in this study are meant 
primarily as a visual aid that provides a bridge between 
the density contour plots and the nearest neighbor plots de- 
scribed below. The basic idea is to connect "segments" of 
occupied cells to give some idea as to the location and ex- 
tent of both connected and island regions in the cluster. An- 
other appropriate terminology might be "tree" structure, but 
as that is used to denote a specific analysis technique, the 
term "segmentation" is used instead. The plot is constructed 
in much the same manner as are the contour plots: cell 
counts of galaxies are tabulated for a specified grid (again, 
0.3x0.3 (/i _1 Mpc) 2 square cells) , the total number of cells is 
increased by subdividing the grid further in order to increase 
the resolution, and new cell counts are computed using the 
bicubic spline interpolation scheme as both an estimator and 
a smoothing mask. The segmentation algorithm is then sim- 
ply the imposition of an additional step that thresholds the 
interpolated cell counts: the count is reset to zero for values 
below the threshold, and a user-defined gray scale is applied 
for counts above the threshold. In practice, we found the 
following choice of thresholds most useful: 

1. For central core densities ao < 100 gal/(/i _1 Mpc) 2 , 
the cutoff is 1.0 galaxy per bin (or, equivalently, 

11.1 gal/(/i" 1 Mpc) 2 ), 

2. For central core densities 100 < ao < 
200 gal/(/i _1 Mpc) 2 , the cutoff is 1.5 galaxies per bin (or, 
equivalently, 16.7 gal/(/i _1 Mpc) 2 ), 

3. For central core densities ao > 200 gal/(/i _1 Mpc) 2 , 
the cutoff is 2.0 galaxies per bin (or, equivalently, 

22.2 gal/(/i _1 Mpc) 2 ). 



The resulting plot exhibits the connectivity of the cluster 
galaxy distribution as a function of cell count threshold. 

A final visualization technique useful for identifying 
clumps that may represent true subclusters is the plot of 
nearest neighbors. The maximum distance used in order to 
define two galaxies as neighbors and the number of neighbors 
required to define a neighborhood of galaxies is, of course, 
arbitrary. In the plots presented in Appendix A, the selected 
distance is 0.25 h _1 Mpc and the minimum number of galax- 
ies is set to four, i.e., each point on the plot is the position of 
a galaxy with four or more galaxies closer to it than 0.25 h _1 
Mpc. Where appropriate, different choices are utilized to 
highlight certain features, and are presented in the detailed 
analyses of individual clusters. 

3.3 The a Test 

The 3-dimensional a test used to det ect substructure i n clus- 
ters of galaxies was i ntroduced bv IWest fc Bothunl Jl990) 
and slightly revised bv lPinknev et alJ (11996ft . The test mea- 
sures the shift in the arithmetic mean centroid of a cluster 
with the galaxy locations weighted by kinematics in spec- 
ified neighborhoods about each galaxy. The centroid used 
is simply the arithmetic mean of all galaxy positions in the 
cluster. The location of each galaxy is weighted by the in- 
verse of the velocity dispersion of nearby galaxies. Thus, a 
neighborhood of galaxies with similar velocities (i.e., low dis- 
persion) has a higher weight and results in a centroid shift 
in their direction. If the member velocities are distributed 
randomly among the member positions, the centoid shift is 
statistically insignificant. 

Thus, the steps leading to construction of the a statistic 
are as follows: 

(i) Convert (RA, Dec) sky coordinates to cartesian co- 
ordinates (x, y), and calculate the centroid of the two- 
dimensional galaxy distribution, 

1 N 1 - 

i=l i=l 

(ii) Assign a weight u>i = l/<Xj to each galaxy, where Oi is 
the line of sight velocity dispersion for galaxy i and its N nn 
nearest neighbors where an arbitrary choice is made to set 
N nn = \/N, 

(iii) For each galaxy i and its N nn nearest neighbors, cal- 
culate the weighted centroid, 

x c = *N nn+1 andy c = '„„„ +1 • (20) 

(iv) The difference in centroid of these nearest neighbor 
velocity groups and the unweighted centroid is calculated, 

0» = \l (Zc - x' c f + (y c - y' c ) 2 (21) 

(v) Finally, the a statistic quantifies the cluster substruc- 
ture as an average of the a.i values, 

1 - 

tt =jvE a - ( 22 ) 
j=i 

The test is normalized by comparison with the value a takes 
for Monte Carlo distributions created by randomly shuffling 
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the velocity among galaxy positions fixed by the data. The 
a statistic for the cluster is compared to the mean a statistic 
from Monte Carlo simulations with the difference expressed 
in standard deviations used to estimate the probability for 
the presence of substructure. A subcluster will influence the 
centroid of its velocity group provided it deviates from the 
global distribution both spatially and in velocity dispersion. 
Detailed results of the a-statistic computed for the clusters 
themselves as wel l as fo r the Monte Carlo simulations are 
contained in lVickl J2004I) . 

3.4 The P Test 

The j3 test is used to detect deviations from mirror symmetry 
about the cluster centre where mirror symmetry is defined 
as a 2-dimensional coordinate inversion (x,y) — > (—x,—y) 
iWest. Qemler. fc D ekcl 1988). After determining the clus- 
ter centroid, the mean distance to each galaxy's nearest 
neighbors is calculated. Then, the location of the inverted 
position of each galaxy is computed, and the mean distance 
from that location to the nearest neighbors is calculated. Ac- 
tual mean distances are compared to the mirror image dis- 
tances. Large differences indicate asymmetry and possible 
substructure. This test is particularly useful in quantifying 
visual impressions of asymmetry conveyed by the visualiza- 
tion plots, and can often be correlated with results from the 
ellipticity algorithm. 

Here, the asymmetry parameter for galaxy i is defined 

as 

A = log (|) (23) 

where do is the mean distance of the galaxy's N nn nearest 
neighbors (N nn = \/~N rounded to the nearest integer), dj is 
the mean distance of the mirror image point of each galaxy 
to the N nn nearest neighbors, and Nrotai is the number of 
total galaxies in the cluster. The /3 statistic for a cluster is 
the mean f3 value for all the galaxies times one thousand or 

P = 1000 (fii) . (24) 

The significance of the P statistic of the data is estimated 
by comparison against a distribution obtained from Monte 
Carlo simulation. This is accomplished by holding constant 
the distance from the centre of each galaxy in the data, but 
uniformly varying the polar angle from to 27r. The galaxy 
coordinates are then relocated at Xi = cos cj> and %ji = 
Ti sin </>. The p statistic is recalculated for each reshuffled 
galaxy cluster. The mean and standard deviation computed 
from the Monte Carlo simulations are then compared to the 
value from the original data. A Monte Carlo probability of 
random occurrence is used to estimate the significance of 
any deviation from mirror symmetry. A significant deviation 
may indicate an asymmetrical distribution caused by one or 
more subclusters. Detailed results for both the /3-statistic 
computed for the clusters as wel l as fo r the Monte Carlo 
simulations are contained in IVickl (l2004Tl . 

3.5 The k Test 

The re test is a three-dimensional substructure test based 
on nearest neighbor velocity relati onships first introduce d 
in a substructure study of A1656 llColless fc Dunnlll99d) . 



As discussed in that study, the re test is similar in intent to 
the Dressler-Shectman A test iDressler fc Shectmanlll988h 
in that it probes deviations in localized velocity variations 
compared to the global cluster velocity distribution. It is 
especially useful for detecting subgroup dynamics in a clus- 
ter core region including the case of a projected subgroup 
centroid overlapping the projected cluster centroid. 

The velocity distribution of a user-selected number of 
nearest neighbors n is compared to the velocity distribution 
of the entire cluster with the test statistic n„ constructed as 

N 

Kn = J2 -k>g[PKs{D > D obs )] , (25) 

i=l 

where the sum is over the galaxies in the cluster. Thus, 
Pks{D > D oos ) is the probability that the K-S statistic D 
is greater than the observed value D bs which is simply a 
standard K-S two sample test with a significance that can 
be computed in a straightforward manner as given in, e.g., 
IPress et al.Nl992h . The re statistic is then the (negative) log- 
likelihood that there is no localized deviation in the velocity 
distribution on the scale of sets of n nearest neighbors. The 
larger the value of re, the greater the likelihood that the lo- 
cal velocity distribution is different from the overall cluster 
velocity distribution. The significance of k„ can also be com- 
puted from Monte Carlo simulations by randomly shuffling 
the velocities of the individual galaxies while holding their 
positions constant. 

For this study, the default number of nearest neighbors 
is set to yN, but the results are relati vely i n sensit ive to 
the precise choice when it is ~ \f~N ; see IVickl (|2004h for a 
detailed comparison between the re and A tests as well as 
the sensitivity of the re test to the selected number of nearest 
neighbors. It is convenient to display the results utilizing 
"bubble skyplots" where at each galaxy position, a circle 
is drawn of radius proportional to the log-likelihood of the 
local and overall distributions not being the same under the 
K-S two-sample test. In this paper, the radii of the re bubbles 
was set to re n /30 for visualization purposes. 



4 DETAILED ANALYSIS OF INDIVIDUAL 
CLUSTERS 

This section provides a detailed analysis of the individual 
clusters from consideration of the results of the statistical 
tests combined with visualization plots. It will be seen that 
while detecting substructure is straightforward, unambigu- 
ous interpretation remains problematic. Further correlation 
with X-ray data is expected to remove some ambiguities, es- 
pecially as relates to cluster membership. The projected (2- 
dimensional) galaxy separations are inferred from the proper 
distance at photon time of emission, d p (t s ). 

4.1 Abell 930 

From Figures fAlf aVfe'l. inspection of Abell 930 surface plots 
shows a low density, moderately elliptical cluster containing 
a (projected) bimodal central region aligned east to west 
with the two components having a projected separation of 
approximately 0.50 h' 1 Mpc. Outside the core there is a 
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Figure 4. Local velocities of the nearest neighbors in Fig- 
ure Ele) for Abell 930. 

relatively high density region (greater than 40 galaxies per 
Mpc) 2 ) located approximately 1.0 Mpc due north 
of the cluster centre. The segmentation and nearest neighbor 
plots, Figures fAlT dWe'). again highlight the two-component 
core region (labels A and B) and concentration of galax- 
ies to the north (label C). As discussed below, these three 
components appear to be 3-dimensional structures as well. 

In spite of the relatively high concentration of galaxies 
to the north, the /3 test yields no indication of global asym- 
metry (Pp — 30%). This is due to the positioning of the 
clustre center between the two components in the central 
regions; positioning the center at the centre of one of these 
components would result in a significant asymmetry signal. 
This re-positioning was not done for this cluster as for others 
in this sample because of the near equality in number den- 
sity and total number of galaxies in each of the components, 
and the lack of an X-ray map that might indicate which of 
the components possesses more mass. 

The ellipticity computation yields an e = 0.29 - 0.45 
over the entire test range from radius R = 0.50 - 1.25 
hT 1 Mpc (Table |2l • As expected, the mean radius of R — 
0.50 ft -1 Mpc results in a position angle aligned approxi- 
mately east to west due to the bimodal core (PA — 75°). 
Average dispersion ellipse radii from R = 0.75 - 1.25 h^ 1 
Mpc are aligned from north to south with an average posi- 
tion angle of PA ~ 10°. We use an average ellipticity and 
position angle e = 0.4, PA — 70° to calculate a core radius 
and central density of R c — 0.47 ± 0.01 and <to = 61 ± 2. 

Inspection of the velocity distribution in Figure IAIIT i 
shows a bimodal profile with the two peaks separated by ap- 
proximately 1500 km/s. Comparison of the members of the 
two spatial core components to the two velocity modes yields 
no relationship, but due to this bimodality, the xtei velocity 
test rejects the hypothesis of normality with > 99.9% CL 
(Table [§. 

The 3D plots and tests also indicate the presence of 
substructure (Figures lAlf fWh'land Table|BJl. In particular, 
the peculiar velocity plot of Figure lAlf g) shows the two- 
component central region in the 2D plots separating into the 



main body of the cluster with peculiar velocities less than 
zero (due to the averaging of the two peaks to arrive at the 
cluster mean velocity) while approximately 15 galaxies have 
peculiar velocities of < 1000 km/s (this is the eastern com- 
ponent seen in the 2D plots). It is easier to distinguish the 
western (A), eastern (B), and northern (C) components in 
the central regions as distinct subgroupings by plotting the 
local velocities of only the nearest neighbors (see Figure 
It should also be noted that the possibility cannot be ex- 
cluded that each of these two larger groupings may contain 
finer-detailed structure not resolved here, or that there is 
additional structure to the north. 

Both the a and k tests indicate some level of substruc- 
ture with probabilities of random occurrence of 3.0% and 
5.8%, respectively. The k bubble skyplot presented in Fig- 
ure lAlf fl shows two areas of possible subclustering: a group 
of seven galaxies show relatively large radii bubbles (indi- 
cating high k statistics) from 0.5 — 1.0 ft -1 Mpc north of the 
centre, and a second group of above average-sized bubbles 
includes galaxies in the western component of the core re- 
gion. Among the seven galaxies to the north with relatively 
high k values, three of the seven are well separated from a 
compact grouping of four galaxies indicated by the arrows 
in Figures lAlf ah (e), and (f) (two have nearly coincident 
angular positions and are indistinguishable in the 2D plots). 
Due to their apparent proximity coupled with the results of 
the k test, and the fact that these four galaxies possess ve- 
locities consistent with cluster membership, we hypothesize 
that they comprise a dynamically bound subgroup belonging 
to the cluster. 

In summary, there appears to be significant substruc- 
ture present in this cluster. The substructure of the central 
region makes Abell 930 a potential site for merger dynam- 
ics, the velocity structure is also potentially significant, and 
there are probably at least one or two bound subclusters 
well removed from the central regions. 

4.2 Abell 957 

Abell 957 possesses a dense core with a core radius R c — 
0.16 ± 0.01 Mpc and a central density ao — 318 ± 1 (see 
Figures lA2f a) , f b~l . (d) ) . In both two and three dimensions, 
the cluster central regions are distinctly bimodal, with a 
second high density concentration of galaxies approximately 
0.5 h~ l Mpc east of the core also possessing a large fraction 
of the total luminosity due to two relatively bright galaxies 
whose velocities are consistent with cluster membership (see 
Figure lA*2T c11. Due to a selection effect in the catalogue, a 
very bright D galaxy residing at the centre of the peak num- 
ber density that is a well-established cluster member (based 
on redshift as well as sky location) was not included in the 
present sample. However, we have added it to the sample, 
and have also re-positioned the cluster center to coincide 
with the location of peak number density. As can be seen by 
comparing Figures fA2f b) and (c), the peak luminosity coin- 
cides with the peak galaxy density at the position of the D 
galaxy. 

With a probability for random occurrence of Pp = 0.5%, 
the /? test detects the asymmetry due to the single structure 
east of the core. The bimodal characteristic gives the clus- 
ter a highly elongated shape (see Tableland Figure^SJa)) 
leading to e = 0.68-0.35 from R = 0.50- 1.25ft" 1 Mpc, and 
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where the semimajor axis is consistently aligned with the 
two central components in this range with a PA — 74° —42° . 
Even at larger distances from the centre, with the outer 
regions being relatively diffuse, the dispersion ellipses are 
strongly influenced by the dense central region. However, 
from Figure [^Jd), it can be seen that the core itself pos- 
sesses some ellipticity with a position angle pointed directly 
toward the eastern subgroup, possibly indicating tidal inter- 
actions between these two central components. 

Considered in its entirety, the velocity distribution is 
consistent with a normal distribution. However, note from 
Figure IA2f a'l that the 0.5 ft -1 Mpc dispersion ellipse re- 
gion contains several galaxies with velocities that are slower 
(open circles) and faster (plus signs) by more than 1.3a com- 
pared to the cluster mean velocity. While the k test for the 
entire cluster shows no significant signal with P K — 14.4%, 
four galaxies in the core show a 5% or less probability that 
the average of their nearest neighbor velocities belong to 
the clusters velocity distribution (see Figures lA2f gWh'O. In 
particular, some of these galaxies may form part of a small 
subcluster with projected (surface) centroid nearly coinci- 
dent with the core centroid (see Figure lA"2T h'Q. and it is the 
relatively slow local velocities giving rise to the large bub- 
bles in Figure lA"2T f), This may be evidence for a line-of-sight 
merger. 

Although the k test does not indicate anything signifi- 
cant in the eastern subgroup, it is not unreasonable to spec- 
ulate that this group may still be moving toward or away 
from the core in a direction perpendicular to the line of 
sight. It should also be noted that Rosat All Sky Survey 
(RASS) data does exist for this cluster, and confirms that 
the cluster's central potential well is very close to the com- 
puted centroid of the core. The region appears dynamically 
active, and also hints that the eastern clump is, in fact, a 
bound subgroup. However, the RASS observations provide 
only 454s of exposure time. This cluster is also at the edge 
(44' off-axis) of a PSPC observation with a total of 1.95 
ksec. A much longer exposure time with better resolution is 
required to infer further dynamical properties. 

In previous st udies of this cluste r, th e bimodal core has 
been i dentified bv lBeers et alJ lll99l j) andjKries sler fc Beers! 
il997t) . The analysis presented in iBeers et alJ (Il99ll) con- 
cludes that the two components are not grabvitationally 
bound, but uncertainty in the radial velocity difference 
does allow for bound solutions. With the larger sample of 
redshifts available here, this question will be investigated 
further in a following paper. Using redshift data for 34 
galaxies, ISolanes. Salvador-Sole, fc Gonzalez-Casadd lll999l'l 
found the velocity distribution to be consistent with a nor- 
mal distribution with statistics similar to those presenterd 
here for 90 galaxies. The Solanes et al. study also performed 
the A test with a result of Pa = 3.4% compared to results 
found here of P A = 17%. 

4.3 Abell 1139 

The analysis of A1139 presented here is based on 106 
galaxies. The fitted optical core radius of Abell 1139 is 
0.26 h~ Mpc with a central density of 138 which is fairly 
typical of the poor clusters in this sample. Figure IA3f b~l 
shows the core and the more asymmetric envelope. The az- 
imuthally averaged fits to the cluster shape are not sensi- 



tive to the relatively low density concentration of galaxies 
seen to the west of the core at (0.75, —0.2) Mpc (see Fig- 
ure IA3f d'l'l. At large radii the ellipticity is low, ~ 0.2, and 
the position angle is ~ 65° . The inner part of this cluster is 
elongated and the ellipticity within 0.75 Mpc is ~ 0.4 with a 
position angle of ~ 100° reflecting the core elongation seen 
in the number density map. It should be noted that there 
are three additional relatively bright galaxies located near 
the centre of this cluster having redshifts consistent with 
cluster membership that are not included in the 2dFGRS 
catalogue. Their inclusion shifts the centre of luminosity to- 
ward the number density centroid compared to the sample of 
106 galaxies used here, but does not alter any of the results 
concerning substructure. 

The two-dimensional /3 test indicates no significant 
asymmetry (Pp = 15%). However, the three-dimensional a 
and k tests indicate the presence of strong substructure, 
each with a probability of < 0.1% that the values are drawn 
from a relaxed distribution. The bubble plot in Figure lXST fl 
shows the two disjoint areas that deviate from the cluster 
mean lying to the west at (1.0, 0.0) Mpc and south (0.0, 
— 1.0) Mpc of the main cluster component. These two sys- 
tems can be seen in the number density plot (Figure^SJb)) 
as low-level extensions from the main cluster. These two 
structures are also seen in the nearest neighbor plots (Fig- 
urc lA3f e'l'l. The segmentation plot (Figure^Jd)) shows the 
core and western subclump clearly while only a hint of the 
southern subclump is evident, but both substructures can 
been seen in the local velocity plot fFigurc fA3f h')'). 

Intriguingly, an inspection of the local velocity plot in- 
dicates that there is a velocity gradient across the cluster. 
This shear seen in the velocity structure may indicate that 
this cluster has a significant rotation. Alternatively, the clus- 
ter may be in the late stages of a double merger with two 
subclumps spiraling in to merge with the main system. 

In contrast to the a and k tests, the velocity distribu- 
tion in Figure lA"3T i) does not show any statistical evidence 
for substructure (See Table 0J. In addition, the velocity dis- 
persion of 484 km/s is somewhat on the low side for the poor 
clusters in this sample. However, the three-dimensional tests 
and the cluster morphology indicate that this cluster is dy- 
namically interesting. This along with the possible presence 
of a velocity gradient in the system indicates that additional 
data from X-ray observations would be useful in clarifying 
the dynamical state. 

Previous investigations of this cluster have 
found little or no evidence for substructure, and 
have not identified the velocity g ra dient found 
here. In particular, [Kriessler^^B eersI il997f) and 
iKrvwult. MacGillivrav. fc Flinl dl999l) concluded that 
the cluster is unimo dal with no substructure while 
IWest fc Bothunl (1990) found marginal evidence for 
substructure. 

4.4 Abell 1238 

With a central density of a = 68 ± 3 gal/ (/i -1 Mpc) 2 , Abell 
1238 is a relatively low density cluster with a somewhat 
irregularly shaped central region of core radius R c = 0.41 ± 
0.02 /i -1 Mpc (see Figures ^Ja), (b), and (d)). In general, 
the cluster has a low ellipticity in the central regions (e = 
0.34 at 0.5 h~ l Mpc) that steadily decreases to insignificance 
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at the largest mean radius of 1.25 ft Mpc. Inspection of 
the luminosity- weighted contours, Figure lA*4T c1. reveals the 
brightest galaxy in this catalogue to the east of the centre 
with cz « 50 km/s different from the cluster mean (21,199 
km/s compared to cz — 21, 148 km/s). 

In general, this cluster shows no statistical evi- 
dence for substructure at greater than the 95% CL. 
These results are consisten t wit h those found by 
iRhee. van Haarlem, fc Katgerd dl99ll) . The projected sur- 
face distribution is relatively asymmetric (Pp = 50%), and 
the velocity distribution is indistinguishable from a normal 
distribution. The large bubbles belonging to the galaxies 
southwest of the core in Figure lA"4T fl most likely do not in- 
dicate significant cluster dynamics. One possible subgroup- 
ing is a small clump of 7-8 galaxies to the east (labeled B 
in Figure lA'lT e')') and slightly south at a distance of about 
1.1 ft -1 Mpc from the centre (see Figures [A4f a) and (d)); this 
grouping also contains the brightest cluster member men- 
tioned previously, and is also tightly grouped in the local 
velocity-position plot of Figure IA4f h'l . A second clump ap- 
pears due south of the centre labeled A in Figure lA"4T e). Both 
of these potential subgroups show some marginally signifi- 
cant substructure as given by the ft test, but a long exposure 
X-ray observation is probably required to ascertain whether 
these are bound subgroups and the extent to which their 
presence affects the overall cluster dynamics. 



4.5 Abell 1620 

At a distance of d p (t e ) = 231 h^ 1 Mpc, Abell 1620 is one of 
the most distant clusters in this study. Inspection of Abell 
1620 2D plots (Figures IA5f al-fell reveals an irregularly- 
shaped, low-density central region with four distinct group- 
ings outside the central region labeled A, C, D, and E 
with a fifth possible subcluster just northwest of the core 
labeled B. Groupings A and C possess maximum num- 
ber densities equivalent to that of the central region, ~ 
60 gal/(/i _1 Mpc) 2 . Groupings D and E are less prominent 
with densities of ~ 30 - 40 gal/(h _1 Mpc) 2 . 

The central regionof the cluster show almost no ellip- 
ticity at a mean radius of 0.5 h -1 Mpc, but are moderately 
elliptical with e = 0.49 at a mean radius of 0.75 ft -1 Mpc as 
grouping B influences the computation. At mean radii of 1.0 
ft -1 Mpc, the ellipticity is still around e = 0.50, but there is a 
position angle change due to the influence of C and E. Over- 
all, the position angles maintain a fairly narrow range at all 
mean radii, ranging from PA — 17° to 41°, in approximate 
alignment with E, the core, and C with the exception of the 
0.5 h' 1 Mpc. The core parameters are computed using the 
average ellipticity and position angle to yield a mean core ra- 
dius of R c = 0.47 ±0.02 and a central density of ao = 59 ±3 
gal/(/i _1 Mpc) 2 indicating the relatively diffuse nature of 
the central region. 

The velocity distribution is skewed toward higher ve- 
locities (Figure fA5f O) resulting in a skewness value of S = 
0.77 = 1.9crs. The distribution is also somewhat broad com- 
pared to a Gaussian with K — 0.65 = 0.6<tk- The %vei test 
rejects the hypothesis of normality at the 97% CL. 

The f3 test for asymmetry indicates a very high prob- 
ability of substructure with a Pp w 0.1%. While the sub- 
groups are somewhat evenly distributed around the centre, 



the density of the subgroups is significantly different, caus- 
ing the computed high asymmetry value. 

The 3D k test indicates a high probability of substruc- 
ture with a P K = 0.9%. As seen in Figure IXST fl. the level of 
K is relatively high in the core, and subgroups C, D, and E. 
The absence of significant k in A can be explained by the 
similarity of A's local velocity distribution relative to the 
global cluster velocity distribution: the mean velocity of A 
is 24,486 km/s compared to 24,430 km/s for the cluster, and 
the related dispersions are 1103 km/s and 1001 km/s. Also 
note the presence of significant k in the (projected) central 
region that may indicate the presence of an additional sub- 
cluster or substructure. 

4.6 Abell 1663 

Surface plots of Abell 1663 (Figures IA6f aWe'0 reveal a 
large, low density, bimodal central region. The centroid cal- 
culation used here is biased by the (projected) bimodal 
structure region so that the centre is re-positioned to be 
coincident with the peak number density. The northeast 
component contains the most luminosity for this particular 
cluster sample. 

Three galaxies are questionable cluster members. The 
second brightest galaxy in the cluster is also the slowest 
galaxy in the cluster (2.9er from the mean) and on the pe- 
riphery of the cluster (1.4 ft -1 Mpc from the centre at co- 
ordinates (—0.72, —1.18)). The other two galaxies are not 
overly bright (bj = 18.6 and 18.7), but are very close to the 
bright galaxy in both space (within 0.19 ft Mpc) and ve- 
locity (within 350 km/s). This proximity suggests a bound 
group of three galaxies that, based upon the spatial, veloc- 
ity, and luminosity data, is probably not associated with this 
cluster. As their influence on the statistical tests is negligi- 
ble, it does not matter if these three galaxies are included 
or removed from the cluster sample. 

Compared to the other clusters in the study, Abell 1663 
is moderately elliptical with a position angle biased by the 
presence of the subcluster to the northeast (see Tabled . The 
core parameter values of R c — 0.46 ± 0.03 and ao = 57 ± 4 
are computed using the ellipticity parameters for the 1.0 hT 1 
Mpc ellipse that improves the fit in this region compared to 
the circular fit while degrading the fit in the outer regions. 

The velocity distribution is skewed toward lower veloc- 
ities resulting in a skewness statistic of S = —0.60 = 1.5crs 
with a kurtosis of K — 0.86 = 0.8o\jf (also see Figure^Ji)). 
However, the xtei test indicates this distribution is indistin- 
guishable from a normal distribution with only a 36% CL 
to reject the normality hypothesis. 

With the re-positioning of the cluster centroid, the f3 
test then yields a marginally significant result of Pp = 4%. 
However, the strong indication from the k test of P K = 0.1% 
together with inspection of the bubbleplot in Figure IXtjT fl 
seems to indicate additional substructure beyond the ap- 
parent bimodality. Two subclusters are evident in the 3D 
plots (Figures |M f M h ))- One appears slightly east of the 
centre and the other approximately 0.8 ft -1 Mpc west of 
the centre. The subcluster nearest the centre is in very close 
proximity to but not actually part of the bimodal central 
region. With the faster subcluster east of the centre and the 
slower subcluster west of the centre, the velocity gradient is 
a possible signature for cluster rotation or infall dynamics, 
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Figure 5. ROSAT PSPC image of Abell 1750 showing the two 
major components are associated with X-ray emission indicative 
of bound clumps. The brightest cluster member is located near 
the peak of the C subcluster emission. 

but further analysis is required for an accurate determina- 
tion. The k bubble plot f Figure lA6f f 1) also shows another 
possible subcluster approximately 1.0 hT 1 Mpc due south of 
the centre. 

4.7 Abell 1750 

With a mean redshift of z = 0.086, Abell 1750 is one of the 
more distant clusters in this study (d p (t e ) = 232 h~ x Mpc), 
and the present sample comprises 78 measured red- 
shifts. It has been investigated for substructure in 
both the g alaxy distribu ti on a n d intracluster X-ray 
emissi o n bv iFo rman et alJ dl98l| ). |Quirrtan^^sRamireJ 
ll98Sh. iRamirez fc Quintanal JlQQOTl. iBeers et al] | l99l [) . 
West. Jones fc Formanl Jl995p , iBliton et alT~ Jl998f) . 
Jones fc Formanl Jl99Sh . and iDonnellv etaD i200lt) . 
Overall, our results are consistent with these previous 
studies although we do find some differences described 
below. 

The cluster shows two prominent galaxy concentrations 
in both the galaxy and the X-ray distributions. In the fol- 
lowing analysis, the cluster center is set to be coincident 
with the peak X-ray emission (labeled as Core in Figure lA7ll 
although the number density in the C subgroup is compara- 
ble (and the visual comparison between the two is affected 
by binning). Fitting any smooth profile over the extent of 
this cluster is problematic due to the clumpy nature of the 
galaxy distribution. For the Core region, the elliptical King 
fit yields the smallest value in this sample for the central den- 
sity of only 41 ± 4 gal/(/i _1 Mpc) 2 . At mean radii greater 
than > 0.50 ft -1 Mpc, the ellipticity calculation is heavily 
biased by galaxies to the north and northeast of the core 
region (see Figures lATT a')- ( d) ) . 

The dumpiness in the projected surface density is espe- 
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Figure 6. Local velocities of the nearest neighbors in Fig- 
ure IA7f e~l for Abell 1750. In this plot, the Core galaxies are 
marked by squares, the A galaxies by diamonds, the B galaxies 
by crosses, and the C galaxies by circles. 

cially evident in Figures fA7f dWe % l where at least 3 potential 
subclusters are indicated by labels A, B, and C with the C 
grouping possibly possessing additional structure. The Core 
a nd C grouping corresp ond to the SW and NE subclump 
of IDonnellv et ail (|^001). The other clumps identified here 
have no corresponding structure in the X-ray maps pre- 
sented by Donnelly et al. or in the ROSAT PSPC X-ray 
map of Figure [£] The peak luminosity for this cluster sam- 
ple is found in subcluster C due to the presence of a single 
bright galaxy approximately 1.0 Mpc northeast of the cen- 
tre. The X-ray emission to the southeast of the core region 
has no counterpart in the galaxy distribution in the current 
catalogue. 

In addition, the velocity histogram in Figure IA7IT l 
shows an additional velocity component corresponding to 
cz ~ 27,000 km/s that is not present in the Donnelly et 
al. data. The hypothesis that the velocity distribution for 
Abell 1750 can be drawn from a normal distribution can be 
rejected at greater than the 98% CL (Table 0}. 

Quantifying the highly asymmetric visual and clumpy 
appearance of the distribution of galaxies on the sky, the j3 
test indicates that substructure is present at the < 0.1%. 
The higher order k test (Figure \A7t ) also shows distinct 
signs of substructure with relatively clear evidence for the 
Core, A, B, and C subclusters (P K < 3%). This test reveals 
that that both the Core and C regions seen in the num- 
ber density map may be composed of 2-3 units. There is 
some evidence for this in the Xray map of Figure for C, 
but the Core shows no obvious substructure in this figure. 
Further support for subclustering is revealed by the plot 
of local velocities of the nearest neighbors shown in Fig- 
ureEJwhere again the Core galaxies may possibly break into 
smaller clumps separated in local velocity by ~ 180 km/s, 
and with clear identifications for the A, B, and C subclus- 
ters. 

Overall, the results above are consistent with previ- 
ous X-ray results that Abel l 1750 is at least a triple sys- 
tem Jjones fc Formanl ll99Sft . The velocity data presented 
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in iBeers et alJ <ll99ll) reveal that this system is even more 
complex with at least four components of which the subclus- 
ters labeled as Core and D here appear to be gravitation- 
ally bound and infallin g. Further dynamical analysis simi- 
lar to that present ed in lBeers. Geller. fc Huchral il982T) and 
IBeers et alJ l)l99lf) applied to the subclusters identified here 
will be prese nted in a future pa per in order to confirm the 
conclusion in lBeers et alJ dl99lT) that the Core and C sub- 
clusters are gravitationally bound, and to ascertain whether 
any of the other subclusters are also bound components. 

4.8 Abell 2734 

Abell 2734 has a well defined, elongated core (see Figures 
EUa) - (e))- With e=0.55, PA= 114° with low to insignifi- 
cant ellipticity at larger radii (Table Comparison of the 
number density contours with the luminosity density con- 
tours, Figures fA8f b~l and (c), shows that the centre of lumi- 
nosity is slightly offset from the number density centroid due 
to the position of the brightest galaxy in the catalogue of 
this cluster (bj = 15.60) at relative coordinates [x — —0.14 
and y = —0.23). Further comparison of the position angle 
of the 0.5 /i -1 Mpc ellipse with the luminosity density con- 
tours in Figure lASf c 1 ) reveals that the line connecting the 
centre of luminosity to the cluster centroid is approximately 
parallel to the position angle of the dispersion ellipse. The 
core parameters are R c = 0.31 ± 0.01 h -1 Mpc with a max- 
imum density go = 121 ± 5 gals/(/i~ 1 Mpc) 2 , and where the 
elliptical fit of the core parameters provides a noticeably 
better representation of the inner regions compared to the 
circularly symmetric fit. 

Inspection of the redshift distribution shows the pres- 
ence of three galaxies with redshifts close to cz = 
22,000fcm/s. Each of those three are possible background 
galaxies, but only one has a magnitude that makes it an 
obvious non-member candidate. These three galaxies have 
an insignificant influence on the computation of the velocity 
statistics that yields relatively high values for the skewness 
and kurtosis of S — 0.71 and K = 1.38. These values are 
2.1cr and 1.6a, respectively, away from the expected value 
for a normal distribution (see Table [1J. However, the xtei 
test yields an insignificant 85% CL that the distribution can 
be rejected as normal. 

Of the four statistical tests, the 2D test yields the 
strongest indication of structure/substructure with the com- 
puted asymmetry having a probability of random occurrence 
of only Pp = 0.4%. Closer examination reveals this asym- 
metry is caused by what appears to be a chance distribution 
- the mirror image of a relatively isolated galaxy happens 
to fall within a small group of galaxies. The member with 
the highest j3 statistic is located near coordinate (0.5, —1.0). 
Its mirror image at coordinate (—0.5, 1.0) is located close 
to a group that also exhibits substructure in the k test (see 
below). Other mirror images with high /3 statistics do not 
fall in groups with any appreciable 3D structure, suggest- 
ing these groups may be the result of projection. If so, the 
positive (3 test result could be a false positive. 

With P a = 44.8% and P K = 35.3%, the 3D tests show 
no significant indication of dynamical substructure. How- 
ever, the k bubble skyplot in Figure lXsT fl shows a small sub- 
cluster in the northeast at approximate coordinates (—0.4, 
1.1), and possibly some structure in the core. The structure 



in the core could be a merger signature, but more detailed 
examination is required to be certain. 

Pri or substructure studies of this cluster found simila r 
results. ISolanes. Salvador-Sole, fc Gonzalez-Casadol lll99STl . 
using 45 members with redshift data, found no substructure 
(Ps = 28.3%, P K = 41.5%, P A = 14.0%). Kolokotronis et 
al. (2001) found weak or no substructure activity comparing 
optical and X-ray s urface brightn e ss par ameters (e, PA, and 
centroid shift, sc) . iBiviano et all {2002), using 77 members 
with redshift data, found no substructure with the A test 
(Pa = 10.0%). The A test result for our sample is P A = 38% 
fsee lVicklEiol) . 

4.9 Abell 2814 

At a distance of d p (t e ) = 284 h' 1 Mpc, Abell 2814 is the 
most distant cluster in this study. Figures ISSfaWb) reveal 
an irregularly-shaped cluster possessing a relatively well- 
defined core region with isolated lower density concentra- 
tions of galaxies well removed from the central regions to 
the east and southeast labelled as groupings A and B in 
Figures lA9f d)- (e) . A few, possibly isolated, galaxies to the 
north and south are also present but do not appear to form 
subgroupings. Comparison of Figures IA9f b'l and (c) show 
that the luminosity density distribution is very similar to 
the number density distribution with most of the light em- 
anating from the central regions. 

The ellipticity and position angle of the 0.5 ft -1 Mpc 
mean radius ellipse are e — 0.30 and PA — 100° with the 
position angle connecting the centres of the regions marked 
CI and C2. Due to the galaxies to the north and south be- 
coming included in the ellipticity computation, the position 
angle shifts by 30° for the 0.75 h' 1 Mpc ellipse (e = 0.27 
and PA = 130°) At larger mean radii, the galaxies to the 
east and southeast again increase the ellipticity while the 
position angle does not change appreciably (see Table In 
fitting the core parameters, the circular fit does well with 
the outer regions while using the ellipse parameters for the 
0.5 ft -1 Mpc ellipse provides a noticeably better fit to the 
inner regions, yielding values of R c = 0.36 ± 0.01 h" 1 Mpc 
and cr = 81 ±4gal/(/i~ 1 Mpc) 2 . 

The velocity distribution shown in Figure IXST il appears 
relatively normal, and the xtei test gives an insignificant 
68% CL to reject the normality hypothesis. Similarly, the 
skew and kurtosis show relatively small deviations from nor- 
mal values (see Table 0J. 

Due to the asymmetry introduced by the concentrations 
to the east and southeast, the (3 test gives a Pp = 8% chance 
of random occurrence. The three-dimensional a and k, tests 
yield P a — 1.3% and P K = 6.5% probabilities of random oc- 
currence, respectively. However, the results from these latter 
two tests are influenced not by the potential substructure in 
the eastern and southeastern concentrations as much as by 
a potential subcluster of at least eight galaxies just west of 
the cluster centre labelled CI (see Figures lA9f f ) and (h)). 
As can be seen from the plot of nearest neighbor local veloc- 
ities in Figure |7[ a), the groups of bubbles in Figure lA"9T f ) 
labeled as CI and C2 appear as a distinct clumps with a few 
remaining galaxies in the projected (2D) core not appear- 
ing to clump together. The brightest cluster galaxy resides 
near the centre of the position-local velocity space. These 
isolated galaxies may reside outside the core region in 3D 
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Figure 7. Local velocities of the nearest neighbors for Abell 2814 in Figure IHT eL For 
visualization purposes, only the galaxies in the (projected 2D) core are shown in (a) with 
the brightest cluster galaxy indicated by the arrow, while only the galaxies in groupings 
A and B are shown in (b). 



space. Whether CI is a subcluster in or near the core it- 
self or is in the foreground remains to be determined, but 
the structure does indicate probable infall/merger dynamics. 
From Figure |3b), the groupings A and B are also probable 
subclusters. 



4.10 Abell 3027/APM 268 

Abell 3027 and APM 268 are essentially the same cluster 
with the primary differences in the 2dFGRS catalogue be- 
ing some galaxies to the east and northeast included in 
A3027 are not included in APM268, and some galaxies to 
the southwest that are included in APM268 but not included 
in A3027 (see Figures lAToT a) and lATTT all. The central re- 
gions are similar but not identical, leading to small but no- 
ticeable differences in the number and luminosity density 
contour plots, Figures [AlOf b). (c) and lAllf b). (c). Anal- 
ysed together, these two clusters present an interesting test 
case for checking the sensitivity of the various algorithms to 
small differences in cluster membership. 

The surface plots for A3027 indicate the central re- 
gion may consist of two components of approximately equal 
galaxy density, but this is not as evident in the plots for 
APM268. With the centre of luminosity residing at the com- 
puted origin in both cases due to the brightest cluster mem- 
ber (6 j = 15.75), no center re-positioning was done (see 
Figures lA"T0T c) and lATiT d). An X -ray map of this cluster 
would be useful in determining the central region distribu- 
tion. The elongated (or bimodal) central structure biases the 
ellipticity computation to yield high ellipticities with posi- 
tion angles oriented north to south for mean radii out to 
~ 0.75 Mpc at which point the galaxies in the outer region 
become included and reduce the calculated ellipticity values 
(see Table |2J. Using the ellipticity and position angle of the 
0.75 Mpc mean radius, the core parameters are computed 
to be R c = 0.29 ± 0.01 and a = 90 ± 3 for A3027 and 
R c = 0.30 ± 0.01 and a = 81 ± 3 for APM268. 

Four groupings of galaxies outside the core are common 
to both clusters, and are marked as A, B, C and E in Fig- 



ures lAlPT d) and lAllF d). A concentration to the northeast in 
Abell 3027 is marked as F in Figure lAToT dl. and a concen- 
tration to the southwest is marked as D in Figure lAl If d). 
Grouping C contains the second brightest galaxy in this 
sample (bj = 16.08) so that the luminosity density is only 
slightly less than in the core. 

Inspection of the velocity histograms in Figures lAlOf i) 
and lAllT i) hint at superpositions of distributions that con- 
sidered as a single entity is clearly not normal. The xtei 
test rejects the hypothesis of a single normal distribution 
at greater than the 99.9% CL and 98% CL for APM268. 
The galaxies in the north contribute most of the slower 
galaxies (marked by the open circles in Figures rXToT a) and 
lAllf a)). and the distribution centered at ~ 23, OOOkm/s re- 
sulting from eliminating galaxies with velocities less than 
~ 22, OOOkm/s is decidely more gaussian-shaped. 

The a and f3 tests yield no significant substructure for 
A3027 but do for APM268 (Table E). With P K = 0.2% and 
2.2%, the k test does show strong indications of substruc- 
ture. The slower galaxies to the north have relatively large 
bubbles, but do not appear to be part of any organized struc- 
ture, and may be in the foreground, and grouping A does 
not appear to be a 3-dimensional structure. Grouping B, just 
outside the core region, does not show significant k, but does 
show some separation from the core in Figures lAlOf f) and 
lAliT f). However, groupings C and E show significant k and 
are possible subclusters. Groupings D and F do not show a 
large amount of k, but appear in the local velocity plots of 
Figures lAlOf h) and lAllT h) (along with C and E), and may 
be subclusters with dominant galaxy motions tranverse to 
the line-of-sight, or subclusters with similar velocity distri- 
butions to the core region. Overall, A3027/APM268 shows 
substructure in the core and outer regions that should be 
further investigated through a detailed X-ray study. 

4.11 Abell 3094 

The 2D plots in Figures IA12I (a)-(e) show a cluster with an 
extended low density central region with two projected sur- 
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Figure 8. Abell 3094 Velocity Distribution of (a) the East Group (x < -0.60 ft -1 Mpc) 
and (b) the Other Cluster Members (x > -0.60 ft -1 Mpc ). 
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Figure 9. Abell 3094 (a) 3D peculiar velocity, and (b) 3D local velocity plots with the 
outlier removed. 



face density peaks oriented north to south. Concentrations 
of galaxies with densities greater than 40 gal/ (ft -1 Mpc) 2 ap- 
pear 0.75 - 1.0 ft -1 Mpc to the east and north-northwest of 
the cluster centre. The luminosity contour plot has a similar 
shape as the high density regions in the isodensity contour 
plot, but without a bimodal central region (at the resolution 
shown here). 

Due to the apparent bimodal central region, ellipticity 
at the smaller mean radii is noticeable with e = 0.55 and 0.49 
at mean radii of 0.5 and 0.75 h~ Mpc, respectively. At mean 
radii of 1.0 to 1.25 ft -1 Mpc, the groupings to the east and 
north-northwest bias the ellipticity computation resulting in 
values of e = 0.37 — 0.33 at postion angles of PA m 140°. 
We use the e and PA for the mean radius of 1.0 ft -1 Mpc 
to compute a mean core radius R c — 0.39 ± 0.01 ft -1 Mpc 
and a central density <to = 88 ± 2 (Table 

The velocity distribution shown in Figure IA12f O is 
skewed toward higher velocities with S = 0.42 (1.15cr from 
normal) while the kurtosis is close to normal. The relatively 
high value of skewness is probably the primary cause that 
the xtei test rejects of the normal hypothesis at a 93% CL 
(Table El . 

The test indicates cluster asymmetry with a 
marginally significant Pp — 10%, and the k test indicates a 
very high probability of substructure with P K < 0.1%. Most 



of the significant k in the core is associated with galaxies in 
the northern projected component. However, while the local 
velocity plot of the nearest neighbors in Figure UHb) does 
indicate some spatial structure in the core, it is not entirely 
clear whether it should be classified as bimodal. Further in- 
spection of the k bubble skyplot in Figure lA12f fl shows that 
most members with a location of x < — 0.60ft -1 Mpc have 
a relatively high k statistic. Analysis shows that this group, 
with the exception of one galaxy, has a significantly slower 
velocity distribution than the other members of the cluster 
with x > -0.60 ft -1 Mpc (see Figure|HJ. 



The mean velocity of the eastern grouping (not includ- 
ing the outlier) is 18,989 km/s or 898 km/s slower than 
cluster members to the west. The velocity dispersion of the 
eastern grouping is a ve i — 401 km/s, significantly smaller 
than the a ve i — 662 km/s dispersion of the galaxies to the 
west (and a ve i = 728 km/s for the entire cluster). For ease 
of visualization, a dashed line is drawn on Figures lA12f aV 
(f) showing the approximate border of the eastern grouping. 
Ten of the 23 members in the east group show a > 5% prob- 
ability that the member and its nearest neighbors are not 
from the same distribution as the entire cluster. When the 
outlier is removed, the number of members with a > 5% 
probability increases from 10 to 18. A possibility for shear- 
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ing or infall motion is also more evident in the 3D plots with 
the outlier removed as shown in Figures [5J a) and (b). 

In the substructure study of ENACS clusters, 
ISolanes. Salvador-Sole, fc Gonzalez-Casadol il999T> used 46 
redshifts to co mpute a A test statist ic of Pa = 0.4%, similar 
to that given in lBiviano et al] (120021) . These res ults a r e con - 
sistent with the value of Pa = 0.4% found bv I Yield ]2004h 
and the k test result reported here for the sample of 108 
2dFGRS redshifts. Other tests in the ENACS study were 
less conclusive. 

In summary, Abell 3094 shows a high probability of 
substructure with the k test, and some indication with the 
Xvei an< l tests (see Table 0. A grouping of 22 galaxies 
located east of x = — 0.60 ft -1 Mpc may represent a fore- 
ground group or a large subcluster with infall or possibly 
shearing motion. 

4.12 Abell 3880 

Inspection of Figures lA"T3T aWb1 shows Abell 3880 as an el- 
liptical cluster with a well-defined, dense core. Comparison 
of Figures lA13f b)-(c) shows that the luminosity density of 
the core region closely traces the number density. From Fig- 
ure lAlSf b 1 ). there are at least three potential subclusters(A, 
B, and C) with two having density greater than 40 galax- 
ies per (ft -1 Mpc) 2 , each of which is just over 1.0 ft -1 Mpc 
from the cluster centroid (see Figure IXilh m. Two of these 
groupings, A and B, show a relatively high level of lumi- 
nosity (Figure lA13f c')'). This combination of high luminos- 
ity with slow peculiar velocity hints at possible foreground 
contamination; inspection of the velocity-magnitude plot re- 
veals at least four galaxies as candidate, but not obvious, 
non- members. 

Closer inspection of the three groupings shows each to 
contain five galaxies in close projected (2D) spatial proxim- 
ity. In particular, grouping B contains the brightest galaxy 
in this sample (b.j = 14.69) with a velocity a little more than 
1.3a away from the cluster mean (one of the non-member 
candidates mentioned above). Of the remaing four galax- 
ies in grouping B, two have velocities very similar to the 
brightest galaxy, « 1050 km/s less than the cluster mean, 
while the other two have velocities as 365 km/s greater than 
the cluster mean. Thus, it is possible that grouping B con- 
sists of a foreground group of three galaxies with projected 
surface positions close to two cluster members. Grouping 
C is similar to B in that of the five galaxies, one is rela- 
tively bright (bj — 15.45) with a velocity more than 1200 
km/s less than the cluster mean, two others have similar ve- 
locities, while the remaining two have velocities more than 
300 km/s greater than the cluster mean. Finally, grouping 
A possesses relatively low luminosity, and consists of three 
galaxies having velocities close to the cluster mean and two 
galaxies having velocities 705 and 917 km/s faster than the 
mean. 

Corresponding to the elongated contours tracing the 
central regions, the ellipticity is relatively high for the 
0.50 and 0.75 ft -1 Mpc ellipses with e = 0.59 and 0.50 
at position angles of P^4 = 150° and 159°. As the mean 
radius is increased to 1.0 ft -1 Mpc, the ellipticity decreases 
to a relatively low value (e = 0.16, PA = 156°) due to the 
influence of the three high density areas outside the core. To 
reduce the influence of the groupings in the determination 



of the core parameters, the ellipticity and position angle 
are averaged for the 0.75 and 1.0 ft -1 Mpc mean radii to 
yield best fit parameters of R c = 0.25 ± 0.01 ft -1 Mpc and 
ao = 161 ± 13 gal/(ft~ 1 Mpc) 2 (a relatively high value for 
the central density). The resulting elliptical fit to the den- 
sity profile is noticeably better than the circular fit in the 
inner regions of the cluster. 

Other than the three groupings described above, there is 
little obvious substructure present in this cluster. The sta- 
tistical a, K, and xtei yield results for substructure below 
the 90 % CL. However, it is worth noting that the near- 
est neighbor plot and the n bubble skyplot indicate that a 
fourth grouping of at least five galaxies, marked D in Fig- 
ure lAlSf e'l. is also a possible subgroup (and might be de - 
tectable in a long exposure X-ray observation) . ISteir] <ll997f) 
initially considered but rejected A3880 for investigation due 
to there being less than 25 galaxies within 0.5 ft -1 Mpc of 
the centre in the catalogue used for that study; in the study 
here, there are ~ 45 — 50 galaxies within a 0.5 ft -1 Mpc circle 
or dispersion ellipse. However, Stein quotes a mean cluster 
velocity and dispersion for the core galaxies consistent with 
what is found here. 

4.13 Abell 4012 

Abell 4012 has a low-density core with one obvious sub- 
grouping of galaxies approximately 1.0 ft -1 Mpc due north 
of the centre with a density exceeding 40 gal/(ft -1 Mpc) 2 
(see Figures fA14f a~l. (b)). These five galaxies in the northern 
grouping are marked as A in Figure IA14I and are discussed 
further below (two galaxies have nearly coincident angular 
positions and are indistinguishable in the plots). 

At small mean radii, the cluster has a moderate level of 
ellipticity that decreases to insignificance at the larger mean 
radii (Table At a mean radius of 0.5 ft" 1 Mpc, e = 0.42 
with PA — 14° where the position angle points directly 
toward Abell 4013. By a mean radius of 1.25 ft -1 Mpc, the 
ellipticity is negligible (e = 0.01, PA = 133°). Compared to a 
circular fit, the elliptical fit is noticeably better with a mean 
core radius of R c = 0.28 ± 0.01ft" 1 Mpc and a = 71 ± 3 
galaxies/(ft _1 Mpc) 2 . 

With a value of a ve i = 471 km/s, the velocity dispersion 
is one of the smallest found in this study (Table 2J. Inspec- 
tion of the velocity distribution in Figure lA14lT ) shows what 
appears to be a somewhat normal distribution. The skew- 
ness and kurtosis values of S — —0.08 and K = 0.37 are 
each less than 0.5<r from what is expected for a normal dis- 
tribution (see Table 2J. However, the Xvei test shows a 96% 
confidence level to reject the normal distribution hypothesis. 

Inspection of Figures lA14f a1-fe1 shows the only clear 
sign of subtructure as the previously mentioned small clump 
of galaxies to the north of centre. However, consisting of 
only five galaxies, this grouping is not enough to cause the 
/3-statistic to indicate significant asymmetry (Pp = 26%). 

The a test indicates substructure [P a = 1% probability 
of random occurrence) but the k test does not (P K = 21%). 
In order to better isolate any apparent dynamical dumpi- 
ness in this cluster, the nearest neighbor positions in Fig- 
urc lAUf c'l are displayed in a local velocity plot in Figure tTHl 
This plot clearly shows the grouping A to be a probable 
bound subcluster. If this interpretation is correct, the lo- 
cal gravitational potential minimum should be detectable 
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Figure 10. Local velocity plot of the nearest neighbor positions 
for Abell 4012 shown in Figure IXUF e'l. 



in a long exposure X-ray observation. The other hints of 
dumpiness in Figure IA14f e) are not as clearly manifested 
in Figure 1101 but cannot be completely excluded as small 
subclusters or merger remnants. 



4.14 Abell 4013 

Inspection of Figures lAlSf al.fb'). and (d) shows a cluster 
with a well-defined, dense core. Approximately 0.5 h^ 1 Mpc 
south of the cluster centre is an additional high density 
group of galaxies having a peak density of greater than 
60 gal/(/i _1 Mpc) 2 . Two small groupings appear west and 
southwest of the centre. For this cluster, the luminosity den- 
sity is very similar to the number density (compare Figures 
IXTBT bl and (c)). 

Within the 0.5 ft -1 Mpc dispersion ellipse, the values 
of e = 0.29 and PA = 163° are biased by the high-density 
group of galaxies to the south. The ellipticity decreases at 
larger radii as the influence of this high-density subcluster 
decreases. At a mean radius of 1.0 h" 1 Mpc, e = 0.20 with 
a PA = 116°. The core fit yilds a core radius of R c = 0.17± 
0.01 and a maximum density of <r = 224±2 gal/(/i _1 Mpc) 2 . 
The resulting fit accurately reproduces the density profile 
for most of the cluster. The f3 test indicates substructure 
at high significance with Pp = 2% where the high level of 
asymmetry is caused by the southern subgroup coupled with 
the relative absence of galaxies to the north. 

Inspection of the velocity histogram (Figure IA15H D 
shows a group of fast galaxies well-separated from the rest 
of the distribution: the majority of the cluster members are 
grouped around a velocity of 15,600 km/s, but a small group 
of 14 are found near 17,600 km/s. This feature explains the 
> 99.9% CL to reject the hypothesis that the distribution is 
normal. These faster galaxies are in three spatial locations, 
and can clearly be identified in Figure IA15f a) (marked by 
'+' symbols). Five galaxies are approximately 1.0 h~ x Mpc 
east of the centre, six galaxies appear near the centre of the 



cluster, and three are approximately 1.0 h~ x Mpc southwest 
of the centre. 

The k test yields an indication of substructure at 
slightly greater than the 95% CL, but inspection of Figure 
IA15f O shows that the large southern subgroup has relatively 
small bubbles, possibly indicating any dynamical interac- 
tion between core and subcluster is perpendicular to the 
line of sight. The high velocity galaxies to the east with the 
large k bubbles are also relatively fainter than the majority 
of cluster members, and may be background galaxies. But 
there are also significant k bubbles in the core region that 
may indicate the merger or dynamical interaction of small 
subgroup (s) (here, spatially unresolved) as indicated by the 
presence of slow and fast galaxies in the core. 

4.15 Abell 4012 and Abell 4013 as a possible 
bound system 

It is worth noting the spatial proximity of these two clusters. 
The projected surface (2-dim) separation between the two 
centroids is » 2.5 h" 1 Mpc, while their mean redshifts are 
cz = 16, 230 ± 498 km/s and cz = 16, 410 ± 904 km/s, re- 
spectively. Figure tTTl shows the galaxy positions and number 
density contours from which it can be seen that the southern 
boundary of A4012 is relatively close to the northern bound- 
ary of A4013. Superimposed on the plot of galaxy positions 
are the 0.5 and 1.0/i -1 Mpc dispersion ellipses showing the 
approximate alignment of the inner regions of the two clus- 
ters. 

Given their proximity, it is worthwhile assessing the 
likelihood that the two clusters are gravitationally bound 
to one another. The cluster masses are estimated using the 
virial relation 

M^gE^feE^) \ (26) 

where N is the number of cluster members, is the pe- 
culiar velocity of the i th member, and ry is the projected 
separation between the i th and j th galaxies. Inspection of 
the velocity distribution for Abell 4013 reveals 14 galaxies 
with velocities significantly faster than the other 71 mem- 
bers, and these are removed in the following analysis in order 
to minimize their effect on the velocity dispersion (see Fig- 
ure 

EUtO). Thus, we estimate the mass of Abell 4012 to 
be « 7 x 10 13 Mq with a velocity dispersion of a ve i = 471 
km/s, and a mass of « 4 x 10 13 Mq for Abell 4013 with a 
velocity dispersion of a ve i — 455 km/s based on eliminating 
the "fast" galaxies. Using the Newtonian binding criterion 

V r 2 Rp < IGMtot sin 2 a cos a , (27) 

and the parametr i c a nalysis described in 
iBeers. Geller. fc Huchral <ll982l) for these masses, pro- 
jection angle a, projected centroid separation, and a radial 
velocity difference of V r = 159 km/s, we find a ~ 95% 
probability that A4012 and A4013 form a bound system 
with a less than 10% probability for the bound 'incoming' 
solutions. The probability calculation is based primarily 
on the possible deviation of the actual radial velocity 
separation being different from the Vr = 159 km/s value 
assuming no worse than a 113 km/s RSS measurement 
error. From Figure 1121 the system is bound as long as 
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Figure 11. Relative galaxy positions and number density contours for Abell 4012 and 
Abell 4013. 
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Figure 12. Bound and unbound solution space for Abell 4012 and 4013. A value of 
Ho = 70 has been assumed for the solution of the parametric equations of motion. The 
angle a is a projection angle between the line-of-sight and the vector joining the two 
cluster centroids. 



V r is accurate to within 200 km/s where the parametric 
solutions indicated in the figu re follow the terminology in 
iBeers. Geller. fc Huchral <1982ft . 



4.16 Abell 4038 

Abell 4038 is a fairly rich cluster of richness class R = 2, and 
we present analysis based on 154 galaxies in this system (the 
most of any cluster in this study). The core of this cluster 
is well defined with a central density of ctq = 261 ± 15 and 



fitted optical core radius of R c = 0.19 Mpc. To the east of 
the core, a dearth of galaxies can be seen in the galaxy dis- 
tribution (Figure |^6l a)) at (-0.5, -0.1) Mpc. This is also 
evident in Figures fA16f b) .(d~l . and (e). This apparent gap in 
the distribution is caused by the bright star (SAO 192167, 
m„=4.6) which limits spectroscopy in the vicinity. The loca- 
tion of the star and a very approximate bounding area inside 
which galaxy identification is probably unobtainable are in- 
dicated by the (circle + cross) in Figures fA16f a~l . (b) .(d) and 
(e). From these figures it can be seen that while the blanking 
area covered by the star represents a substantial fraction of 
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the gap area, there are still portions to the north and south 
that are clearly devoid of galaxies. If the gap is not real, 
one expects to see galaxies just outside the blanking area, 
so that it appears that at least part of the galaxy-free region 
is a true feature of the distribution. But it should be noted 
that the inability to detect cluster members in the vicinity of 
SAO 192167 could significantly bias the results of any sub- 
structure tests using the global cluster galaxy distribution 
on the sky if, in fact, cluster members are not being included. 
However, substructure tests that evaluate local parameters 
or the velocity distribution should be less affected. 

Not surprisingly, the gap in the galaxy distribution 
strongly affects the analysis of this cluster including the el- 
lipticity computation. Looking at the smallest radius of 0.5 
h _1 Mpc, the void has a strong influence on the ellipticity 
by cutting off the distribution on the east side. This results 
in a relatively large ellipticity of 0.43 and a position angle 
of 21°. At the largest radius the ellipticity has risen to 0.50 
and the postion angle has moved to 96°. This is due to the 
apparently disjoint galaxy clump at ( — 1.5, 0) Mpc and the 
transition would possibly be more gradual without the data 
gap. 

The a and (3 tests indicate the presence of strong sub- 
structure. Since these two tests are more heavily influenced 
by the global distribution of galaxies on the sky, the gap in 
the galaxy distribution biases the results of these tests. The 
k test also indicates the presence of substructure with a very 
high significance but since this is a local test, it should be 
less influenced by the gap in the data. However, the abso- 
lute significance of the test needs to be investigated in more 
detail for this cluster. 

The bubble plot in Figure IA16f f) shows at least four 
areas with evidence of substructure. The most obvious are 
two groupings ~ — 1.2 Mpc to the east, and separated from 
the main cluster by the gap discussed above. However, these 
features show up in the velocity plots (Figures IA16f g'l and 
(h)), and are clearly at discrepant velocities compared to 
the main cluster. The next most discrepant grouping is to 
the south of the main cluster at (—0.4, —1.0) Mpc, but this 
grouping is only weakly apparent in the number or luminos- 
ity density plots. This subclump appears to be well sepa- 
rated from the main concentration, but may be affected by 
the gap near the bright star. One other subclump is located 
west of the centre. The clump at (1.0,0.0) Mpc shows up as 
a low-level contour in the number density map and also can 
be seen in the local velocity plot. There may be additional 
substructure just northeast of the core (see Figures fA16f b')') 
and (d)). 

Although the a, /3, and n tests yield strong signals for 
substructure, the xtei test rejects the normality hypothesis 
at only the 96% CL. Yet the velocity distribution as shown 
in Figure lAl6f i) does show hints of being a superposition of 
two (or more) distributions. 

This cluster has been previ ously studied for su b struc- 
ture in the galaxy distribution bvlGirardi et alJ <ll997f ). ISteinl 
Jl997l) . and lKriessler fc Beers ! Jl997l) . and using ROSAT X- 
ray data bv lS^nuecker^t^l .1 (1200 ll) . However, these studies 
considered the central region of the cluster, and did not in- 
clude the eastern extension that may or may not be associ- 
ated. The velocity statistics and A test value are consistent 
with those found here for that region. The Kriessler & Beers 
study presents a contour plot of the central region signifi- 
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Figure 13. Nearest neighbor local velocities versus the y- 
coordinate for Abell S141. The apparent clumping in 2D also 
appears in 3D for at least groupings A, B, and C. 



cantly different in appearance from ours, and their KMM 
analysis identifies three modes that we are unable to re- 
solve. However, a new study of this cluster in X-ray with 
longer exposure times should further clarify the true sub- 
structure and resolve the problematic constraints imposed 
by SAO 192167. 



4.17 Abell S141 

At a distance of d p (t e ) = 57 h' 1 Mpc, Abell S141 is the 
nearest cluster in the study (cz = 5793 km/s). Inspection 
of Abell S141 2D plots in Figure lA"l7l shows a highly elon- 
gated structure with two apparent dumpings in the south 
and possibly one in the north. The core is relatively dense, 
but, together with S1043, S141 has the lowest luminosity of 
the clusters studied here. 

For the four mean radii from 0.5 to 1.25 h~ x Mpc, the 
ellipticity and position angle have a narrow range from e = 
0.67 - 0.58 and PA = 0° (orl80°) - 168° (see Table 0. We 
use the average ellipticity of e = 0.64 and position angle of 
PA = 173° to calculate a core radius R c = 0.26 ± 0.02 and 
a central density 00 = 158 ± 18. 

The velocity dispersion of a ve i = 403 km/s for the dis- 
tribution shown in Figure IA17ID is the lowest in this study. 
Also from this figure, it is evident that the distribution is 
skewed toward higher velocities due to the ten relatively fast 
members labelled with '+' symbols in Figure lATrT al. The 
Xvei rejects the normality hypothesis at the 98% CL, but it 
is difficult to correlate the results of this test as due to any 
of the groupings except for possibly the one farthest south. 

The three apparent subclumpings in the segmentation 
and nearest neighbor plots are labelled as A, B, and C in 
Figure lA17f e~). The local velocities of the nearest neighbors 
shown in Figure [T3l appear to support the hypothesis that 
groupings B and C form subclusters with possibly also a 
fraction of grouping A forming a subcluster. However, the 
statistical tests yield no clear signal for substructure with 
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the p and a tests yielding P = 23% and P a = 34.6%. The 
lack of a signal in the test is easily explained from the ob- 
vious symmetry in the projected directions. The low velocity 
dispersion combined with the 2D symmetry is probably re- 
sponsible for the lack of a signal from the a test. 

Although the k test result of P K = 19.8% using ten 
nearest neighbor velocities (~ v~N) is not statistically sig- 
nificant, there are some indications of subclustering in the 
bubble skyplot in Figure lAl7f f) corresponding to the group- 
ings described above. Performing the n test with numbers of 
nearest neighbor velocities other than y/N did not apprecia- 
bly change the significance level. However, for this cluster, 
the small dispersion of velocities naturally reduces the dis- 
criminating power of the test. This is probably due to the 
geometrical effects introduced from analyzing a cluster with 
one relatively long and two short symmetry axes along a 
direction perpendicular to the long axis. Thus, we conclude 
that in addition to the possible subclustering, the global 
shape of this cluster is triaxial with one long and two short 
axes (or possibly prolate). This geometry may be linked with 
the supercluster environment in which S141 formed. 

4.18 Abell S258 

Abell S258 is highly elongated, as is evident upon viewing 
the 2D plots (Figures IA18f a~) - (e)). Its length stretches over 
3 ft -1 Mpc in an east to west direction while its width is 
typically 1 ft -1 Mpc. It has a moderately dense central re- 
gion that also appears highly elliptical in an east to west 
direction. 

The ellipticity and position angle are very consistent 
for all mean radii with ranges of e = 0.75 - 0.58 and PA 
= 114° - 104° (see Table 0. This cluster has the highest 
ellipticity in the study at a mean radius of 0.5 ft -1 Mpc and 
the highest or second highest ellipticity at mean radii from 
0.75 - 1.25 ft -1 Mpc. Using the ellipticity and position angle 
for the 1.0 ft -1 Mpc test radius (e = 0.71, PA = 104°), we 
find a mean core radius, R c = 0.28 ± 0.04, and a maximum 
density, a = 108 ± 1 (Table 0. 

The velocity distribution appears bimodal (see Figure 
IA18IT )'). This characteristic, however, is not sufficient to 
cause a negative test of normality as the Xvei test rejects 
the normal hypothesis at an insignificant 70% CL (Tablc^J. 

Only the k test indicates significant substructure (Table 
|EJ, but that test results in a very high significance level with 
a P K < 0.1%. Inspection of the k bubble skyplot (Figure 
IA18f f1') shows a large number of members with a high level 
of k in the far eastern part of the cluster. Of the 14 members 
east of x = —1.0 ft -1 Mpc, 12 are grouped very close in 
position and velocity. These 12 galaxies are found to have 
a much slower mean velocity (v = 16,863 km/s), are more 
narrowly dispersed (<j ve i = 272 km/s) compared to the other 
members of the cluster (v = 17,584 km/s and a ve i = 526 
km/s), and are responsible for the bimodal appearance of 
the velocity distribution. They can be seen clearly in the 3D 
velocity plots, Figures lAlSf g) and (h). 

It is unclear if this eastern sublcluster is, in fact, part 
of Abell S258. It is possible that this group is a small fore- 
ground cluster. The eastern group has an average redshift 
of 767 km/s less than the remainder of the cluster which 
could equate to a foreground distance of 11 Mpc (using h 
— .70). However, the velocity difference of 721 km/s is also 



Table 5. Abell S258 Comparison of East and West Group. 



East Group West Group 



Cluster 



N 

V 
a ve l 

bj 



12 

16,863 km/s 
272 km/s 
18.0 



75 

17,584 km/s 
526 km/s 

18.0 
CL = 67% 



87 

17,485 km/s 
557 km/s 

18.0 
CL = 70% 



within the range of possible orbital velocities estimated for 
a subcluster orbiting the main body. A binding probabil- 
ity analysis will be presented in a future paper. Comparing 
velocity versus magnitude for each member in each group 
is inconclusive in determining spatial separation. Moreover, 
removing the eastern group members did not significantly 
improve the likelihood for normality of the velocity distri- 
bution. Table compares key statistics of each group and 
the cluster. 

There are also possible subclusters to the west and 
north that are more readily apparent from inspection of the 
nearest neighbor plot, Figure lA18f e'l. in which they are la- 
beled 'A' and 'B'. Subcluster B appears to contain eight 
galaxies grouped in both position and velocity space with 
other galaxies in the surface plot probably not being sub- 
cluster members. The northern A grouping appears to be a 
result of chance projection in the surface plots. Finally, the 
segmentation, nearest neighbor, and 3D local velocity plots 
suggest the presence of a possible subcluster just southeast 
of the core that also contains the brightest cluster member. 
It is interesting to note that all of the apparently real group- 
ings/subclusters lie (roughly) along the position angle of the 
major axis of the dispersion ellipse. 



4.19 Abell S301 

Abell S301 (DC-0247-31) is one of the nearest clusters in the 
study with an estimated distance of d p (t e ) = 65ft _1 Mpc 
(cz = 6669 km/s). The surface plots in Figures lA19f aWe) 
show the central regions to contain a compact high density 
core and a smaller second component to the southeast. Com- 
parison of the luminosity and number densities in the central 
region shows them with similar structure (Figures IA19f b~l 
and (c)). In addition, one other subgrouping having density 
greater than 40 gal/(/i _1 Mpc) 2 appears to the southwest, 
approximately 1.2 h~ Mpc from the centre, and a lower 
density grouping appears to the northeast at approximately 
1.0ft -1 Mpc. 

Biased by the southeast component, the ellipticity at a 
mean radius of 0.5 Mpc is e = 0.36 at a position angle of 
PA — 161°. As the mean radius increases, the southwestern 
and northeastern groupings influence the calculation to shift 
the position angle with, for example, e = 0.61 and PA = 
48° at a mean radius of 1.0 Mpc (Table |5J. The variation 
in position angle is problematic in achieving an elliptical 
core fit accurate at all radii, but using the average ellipticity 
parameters from the 0.75-1.0 Mpc mean radii is somewhat 
better than using a circular fit and yields core parameters 
of R c = 0.24 ft" 1 Mpc and a = 133gal/(ft _1 Mpc) 2 . 

Inspection of Figure IA19IT I reveals a distribution with 
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most of the galaxies broadly distributed between ~ 6300 — 
7000 km/s and with relatively small tails resulting in a kur- 
tosis of K — 1.9 = 2.0ai{. The Xvei t es t rejects this as a 
gaussian distribution at the 99.9% CL. Two of the slowest 
galaxies (cz « 5000 km/s) are relatively bright compared 
to the other cluster members and are possibly in the fore- 
ground. 

While the a test gives no indication for substructure 
(Pa = 24.2%), both the /3 and k, tests yield statistically sig- 
nificant signals (Pp = 2% and P K = 0.4%). Inspection of 
the bubble and local velocity plots (Figures IA19f f) and (h)) 
reveals a subcluster in the main component of the core not 
evident in the surface plots. In particular, Figure IA19f h) 
shows that the local velocities of these galaxies are signifi- 
cantly faster compared to the cluster mean velocity. Overall, 
the core galaxies have an average velocity of 6968 km/s com- 
pared to the mean value for the cluster of 6594 km/s, but 
13 galaxies in the core have local velocities exceeding 7200 
km/s. These galaxies reside in a very close group near the 
centre of the main component (x = —0.10 to 0.10 h -1 Mpc, 
y — —0.02 to 0.11 hT 1 Mpc). The northeastern grouping also 
appears in the bubble skyplot, but the farthest southwestern 
grouping is virtually absent. Assessing whether the central 
velocity structure is truly in the core or is a result of pro- 
jection effects is required to further determine this cluster's 
dynamical configuration and state. 

This cluster has been investigated previously for 
substructure several times, but with significantly fewer 
redshifts than the 95 us ed here. Using 26 redshifts, 
iDressler fc Shectmanl lll988tl found no sub s tructu re (Pa = 
34.5%). Using 29 redshifts. lEscalera et all (Il994l) found no 
substructure from app lying the A a nd Lee tests (Pa = 52.2% 
and P Lee = 80.0%). ISteinl l|l997l ) performed the A, Lee, 
and velocity normality tests using 25 members and found 
no su bstructure (Pa = 29%, PLee = 63%, and Pi sr orm = 
43% ). |palvador-Sole. Sanroma. fc Oonzalez-Cas adol dl993T ) 
and ISalvador-Sole. Gonzalez-Casado. fc Sola nes ( 1993|) us- 
ing autocorrelation tec hniques found signs of su bstructure. 
Using a KMM method, iKriessler fc Beersl (|l997T l found the 
central region to be bimodal (P = < 0.1%) with the same 
mode centre locations identified here. 

4.20 Abell S333 

Abell S333 is one of the most diffuse clusters in the sample 
with a low surface density core region (see Figures lA20f a) 
and IA20f b'l'l. The luminosity- weighted contour plot (Fig- 
ure [^2lc)) shows the brightest area in the cluster to be 
separated from the core by approximately 1.0 /i _1 Mpc (pro- 
jected) northwest of the centre. The high luminosity is due 
to the brightest galaxy in the cluster (bj = 15.4), but with 
its velocity being several hundred km/s faster than the clus- 
ter mean, it seems unlikely this galaxy is in the foreground. 

The cluster shows consistent ellipticities and position 
angles at all radii tested (see Table 0. The ellipticity ranges 
from e = 0.59 - 0.41, and the position angle range of PA = 
131° - 156° is consistently aligned with the high luminosity 
region (Figures IA20f a~l and (c)). Given the consistent e and 
PA at all distances from the centre, we use the average values 
(e = 0.5 and PA = 135°) to calculate the core parameters. 
The resulting values are R c — 0.35 ± 0.01 and oo = 62 ± 2. 

Inspection of the velocity distribution in Figure |A20IT | 




Figure 14. 3D Local Velocity Plot for Abell S333. This is Figure 
IA20f h'l rotated counter-clockwise 50° about the z-axis. This view 
clearly shows the spatial asymmetry in the local velocity with the 
northern half of the cluster (the leftmost two quadrants in the 
figure) moving away from the line of sight and with the southern 
half (the rightmost two quadrants) moving toward the line of 
sight. 



shows a group of relatively slow galaxies in addition to the 
unusually-distributed main group. The slow group consists 
of eight galaxies that have velocities grouped near 17,500 
km/s compared to the cluster mean of 19,372 km/s with 
their positions being indicated by the 'o' symbol in Figure 
lAloT al. Although moving significantly slower than the clus- 
ter mean velocity, these galaxies are not obvious foreground 
contamination: six have magnitudes greater than the cluster 
mean of 18.1. Two of the galaxies are both brighter than av- 
erage (bj = 16.9 and 17.8) and significantly slower than aver- 
age (v = 17,361 and 17,277 km/s). While these two galaxies 
are possible non-members, their local motions are consistent 
with neighboring galaxies that are members (see below). 

The histogram of the velocities of the majority of the 
cluster members shows them to be unusually distributed 
about a narrow peak near the mean velocity. This is re- 
flected by skewness and kurtosis values that are 2.5cr and 
0.9a away from normal values, respectively (see Table 
The Xvei velocity test rejects the normal hypothesis with a 
99.9% CL. If the eight slowest galaxies are removed from 
the computation, the S and K statistics move closer to nor- 
mal; however, the xtei velocity test still rejects the normal 
hyposthesis at a 93% CL. 

The 3D plots and tests show signs of substructure with 
the k test having the strongest indication with a probability 
of random occcurence of P K < 0.1% (Table |SJ. The sub- 
structure detected by the k test is found in the high lumi- 
nosity region to the northwest of the core as seen in the k 
bubble skyplot (Fi gure lAltjlm The 3D plots also highlight 
this area (Figures IA20f g) and (h)) where ten galaxies in 
the northwest quadrant show local velocities above the 95% 
level of the K-S statistic (represented by the filled circles). 
The spatial asymmetry of the local velocity plot is seen more 
clearly when Figure lA20f Ii) is rotated counter-clockwise 50° 
about the z-axis as shown in Figure HTI 
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(a) (b) (c) 




Figure 15. Galaxy positions for Abell S1043 showing the locations of (a) all galaxies, 
(b) galaxies with v < ll,000km/s, and (c) galaxies with v > ll,000km/s. 
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Figure 16. Abell S1043 K bubble skyplot with v > 11,000 km/s 
members removed. 



To summarize, Abell S333 has a consistently elongated 
shape with two relatively luminous regions (the core and 
the region to the northwest). The luminous region outside 
the core is also where most of the velocity substructure is 
found. Together with the (approximate) north-south asym- 
metry in the local velocity, the cluster shows evidence for 
dynamics possibly involving rotation, shear, or infall within 
its supercluster environment. 



4.21 Abell S1043 

Inspection of the scatter and number density contour plots 
of Abell S1043 shows a cluster with a well-defined, relatively 
dense, elongated central region (Figures IA21l (a'l and (b)). A 
large number of galaxies at distances greater than 1.0 h~ x 
Mpc from the centroid are located to the south and south- 
west of the cluster central regions. As most galaxies to the 
north are less than 0.5 h~ Mpc from the east-west centre 
line, this gives the cluster an appearance of asymmetry. The 
luminosity contour plot (Figure IA21I fell shows two areas 



outside the core with luminosity the same or greater than 
the core. The two areas of greatest luminosity appear in the 
galaxies to the south and southwest mentioned above where 
the four brightest galaxies are located, but with velocities 
showing no evidence for obvious non-membership. 

At mean radii of 0.5 to 0.75 Mpc, the ellipticity 
ranges from e — 0.60 to 0.43. with a position angle from 
PA = 152° to 169° (Table^J. At a radius of 1.0 h' 1 Mpc, 
the galaxies to the east and west of the core tend to reduce 
the ellipticity (e = 0.03) with no shift to the position angle 
(PA = 162°). When the mean radius is increased to 1.25 
h~ x Mpc, the galaxies to the south and southwest influence 
the calculations such that the ellipticity increases to e = 0.32 
and the position angle shifts to PA — 16° . To fit the central 
regions as accurately as possible, the ellipticity and position 
angles for the 0.5 and 0.75 hT 1 Mpc dispersion ellipses are 
used to find a core radius of R c — 0.31±0.01 and a maximum 
density of oo = 107 ± 3. 

From Figure |A2 If i), the velocity distribution is highly 
skewed (the value of S = 0.86 is 2.3a away from the ex- 
pected value of zero) with a "slow" group of galaxies between 
9,000 and 11,000 km/s having a normal-shaped distribution, 
and with an almost uniform distribution of "fast" galaxies 
from 11,000 to 13,500 km/s. For the entire distribution, the 
Xvei test rejects the normal hypothesis with a CL > 99.9%. 
However, when the fast galaxies are removed, the remaining 
galaxies have a velocity distribution indistinguishable from 
a gaussian (CL = 29% to reject). Further analysis reveals 
that relatively few members of this fast group are located 
within the 0.5 h~ x Mpc mean radius ellipse (see Figure rToll . 

The two-dimensional plots and tests are dominated by 
the galaxies to the south and southwest as shown by the seg- 
mentation plot (Figure lA21f d1'h the nearest neighbor plot 
(Figure Elit e)), and the (3 test result of Pp = 0.7%. Com- 
parison of Figure [TBI with Figures [A21f b1. (d), and (e) show 
that ten of the faster galaxies are located nearly due west 
of the central regions, but have only marginal density or 
connectedness relative to the other concentrations. 

Before partitioning the galaxies into the two veloc- 
ity groups, the 3D tests plots indicate a marginal to high 
level of substructure: P a = 10.2% and P K < 0.1% (Ta- 
ble |SJ. However, when the fast group of galaxies mentioned 
above is removed, the significance is substantially reduced 
to P a K low = 8.9%. Comparison of the k bubble skyplot for 
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Figure 17. The central region of APM 917 in 2D and 3D Space. 
Plot (a) shows the central region in 2D space while plot (b) shows 
the central region in 3D as viewed from the west, using local veloc- 
ities for the z direction. Members of component A are denoted by 
squares, component B by circles, and component C by diamonds. 



are no obvious subgroupings in the surface distribution, it 
is expected that the results of this test are due to the pres- 
ence of significant velocity substructure in the core. A more 
detailed examination of the central region in the 3D surface- 
velocity space supports this conjecture. 

The central region appears to be composed of three 
components as seen in Figures Il7f a~l and (b). We classify 
galaxies into one of the three components (A, B, or C) based 
on their 3D proximity to the component. A certain amount 
of arbitrary judgment is used, and the component bound- 
aries are not meant to be necessarily absolute. When viewed 
from the west, it is clear that the A and C components are 
moving slower than the central component, B. All three com- 
ponents have a rather small dispersion among themselves. 
A random reshuffling of member velocities would result in 
higher dispersions resulting in lower weights for the compo- 
nent members and, in turn, a higher centroid shift in the a 
test. 



the slower group in Figure fTTTl with Figures lA21f fl. (g), and 
(h) reveals that it is the fast galaxies to the west that heav- 
ily influence the k test results for all the galaxies. However, 
even with the fast galaxies removed, the central regions still 
show indications of substructure dynamics. In particular, 
there appears to be a possible subcluster near the core that 
is the source of the relatively large remaining k bubbles in 
Figure [Pol 

It is interesting to note that the dispersed galaxies to the 
south and southwest are not inconsistent wi th post-merger 
asymm etries found in the simulations of IPinknev et al.l 
( 1996). However other explanations are possible, and a more 
detailed analysis of the cluster is required. 

4.22 APM 917 

APM 917 appears to have a high-density, well-defined core 
(Figures IA22f a~l - (e)), and does not appear to have regions 
of high-density galaxy groupings outside the central region 
with the luminosity being strongly concentrated in the cen- 
tral region (see Figure lA22f c')'). The velocity distribution ap- 
pears normal. This appearance of normality is further sup- 
ported by the skewness and kurtosis statistics (S = 0.18 and 
K — 0.03) as well as the \vel test yielding only a 73% CL 
to reject the normal hypothesis (see Table 2J. 

The central region, while well-defined, has an unusual 
triangular or almost crescent shape. It has a relatively mod- 
erate level of ellipticity (Table |5J at all mean radii with 
varying position angles (e = 0.45 - 0.27 and PA — 17° - 
56° for R = 0.5 - 1.25 h' 1 Mpc). The average of the el- 
lipticity and position angles for the 0.75 and 1.0 h -1 Mpc 
mean radii are used to calculate the King profile core radius 
and maximum density (Table [3J. The cluster core radius at 
R c = 0.14 ± 0.01 ft -1 Mpc is the smallest in the sample and 
the maximum density, <ro = 270 ± 1, is second highest. 

The results of the k test suggest a possibility of sub- 
structure with P K = 9.0% (Table [fjj. Inspection of the k 
bubble skyplot in Figure lA22f fl shows four members to the 
east of the centre with a high level of ft with some lower 
k values in the central region. The a test indicates a high 
probability of substructure with P a = 0.6%. Since the a 
test measures the velocity dispersion-weighted centroid shift 
when the cluster velocities are randomly shuffled, and there 



4.23 APM 933 

Inspection of APM 933 galaxy position and number density 
contour plots, Figures fA23f a~l and (b), shows a high-density, 
somewhat elongated core, the centre of which is offset from 
the cluster mean centroid by ~ 0.1 Mpc. A grouping of 
galaxies with a density greater than 40 galaxies/ (ft -1 Mpc) 2 
is located approximately 1.3 ft -1 Mpc to the north- northeast 
of the centre. This grouping is also readily apparent in the 
segmentation and nearest neighbors plots, Figures IA23f d) 
and (e). Comparison of the number density and luminosity 
density contour plots, Figures |A23f bland (c), shows that 
the light and mass do not exactly overlap for this cluster. 
Including the grouping to the northeast mentioned above, 
there are three areas to the north with luminosities equal- 
ing or exceeding that of the cluster central region. However, 
correlation of the velocity-magnitude plot with these high lu- 
minosity regions reveals no obvious foreground candidates. 
In fact, the brightest galaxy in this cluster catalogue (at 
the position farthest west in Figure lA23f 0*0 has a measured 
velocity of more than 1.0<r faster than the cluster mean in- 
dicating that it is either an exceptionally bright background 
galaxy or possesses unusual dynamics. 

In general, the ellipticity algorithm converges around 
moderate to low ellipticity solutions. The elongated appear- 
ance of the core is quantified by the ellipticity computation 
for the 0.5 and 0.75 h' 1 Mpc mean radius ellipses yielding 
yielding values of e = 0.53 - 0.40. The fitted core parameters 
are R c = 0.36 ± 0.01 and a = 91 ± 2. 

From Figure |A23IT |. the velocity distribution contains 
a peak at « 14, 600km/s and an extended tail of galaxies 
with velocities greater than 16,000 km/s. These two features 
result in skewness and kurtosis values of S = 0.77 and K — 
0.81 that are 2.2crs and 0.9o\r- , respectively, away from the 
values expected for a normal distribution. The skewness and 
kurtosis combine such that the Xvel test rejects the normal 
hypothesis with a CL of 97%. 

With a value of Pp < 0.1%, the /3 test indicates a high 
probability for nonrandom asymmetry. Although the a test 
does not indicate a high probability for nonrandom sub- 
structure (P a = 18.5%), both the k test (P K = 2.6%) and 
the local velocity plot (Figure rA23f h')') reveal several poten- 
tially interesting subgroupings and possible bound subclus- 
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ters labeled A through D in Figure IX23f e1. The 3D plots 
again show at least eight of the galaxies in the grouping 
to the north-northeast identified as D probably reside in a 
bound subcluster (Figures IA23f fl-(h')'l. In the 3D local ve- 
locity plot, seven of the galaxies in this subcluster have a < 
5% or less probability that the velocity distribution of their 
nearest neighbors belong to the cluster velocity distribution. 
The k bubble skyplot and the local velocity plot also show 
substructure in the central regions highlighted by grouping 
B that is suggestive of two or three subclusters having mo- 
tions consistent with merger dynamics or infall, but a more 
detailed analysis is required for unambiguous identification. 
A smaller group of six galaxies is found outside the central 
regions, and is denoted in the nearest neigbors plot as C. 
Four of these galaxies have velocities within a narrow range 
of 500 km/s, and likely comprise a small subcluster while the 
other two galaxies probably do not belong to this grouping. 
Grouping A does not appear to be a dynamically significant 
subcluster, and may be a binary with one or two additional 
galaxies seen in projection in the 2D plots. 




Y (Mpc) 



Figure 18. EDCC442 local velocities of the nearest neighbors 
shown in Figure lA25f e~). Groupings B, C, D, and E are possible 
subclusters, and there is velocity structure in the core region. 



4.24 EDCC 365 

Inspection of EDCC 365 scatter, isodensity, segmentation, 
and nearest neighbor plots (Figures IA24f a1. (b), (d), and 
(e)) shows a relatively diffuse cluster with a well-defined 
core, a large clump of galaxies just northwest of the core 
(A), and some structure south and southeast of the core that 
does not appear to be subclustering. The luminosity contour 
plot, Figure lA24f ci. highlights the luminosity of the core al- 
though the dominant luminosity signal is in the south due 
to the presence of the brightest cluster galaxy (bj = 15.1). 
The velocity of the brightest galaxy is consistent with clus- 
ter membership (cz = 17, 868 km/s compared to the clus- 
ter mean redshift of cz = 17, 773 km/s). The approximate 
north-south alignment of the groupings biases the ellipticity 
algorithm to yield relatively high ellipticities of e ~ 0.5 — 0.7 
at position angles PA ~ 150 — 160° (see Table |2J. We com- 
pute the core parameters by averaging the ellipticities and 
position angles from the 0.5 and 0.75 h -1 Mpc results to find 
a core radius and maximum density of R c = 0.37 ± 0.01 ft. -1 
Mpc and a = 70 ± 2. 

The velocity histogram fFigurc lAUm appears normal. 
The skew and kurtotsis do show some signs of anormality 
(S — 0.72 and K — 1.00) however, the xtei test has only a 
83% confidence level to reject the normal hypothesis. 

The cluster asymmetry measured by the /3 test yields a 
marginally significant Pp — 4% probability of random occur- 
rence, while the a test gives no indication of substructure 
(P a = 39.1%). However, with P K = 0.2%, the k test in- 
dicates the presence of substructure at a high significance 
level. From the bubble plot of Figure IA24f fl , the substruc- 
ture appears in two distinct regions corresponding to the 
core and the galaxies to the southeast. These galaxies are 
also readily apparent in the peculiar and local velocity plots 
in Figures IA24f g) and (h) , but do not form a well-defined 
clump. Correlation of the galaxies in the core with the larger 
bubbles and the x-y-local velocity plot reveals no obvious 
subclustering in this region. The velocities of the galaxies in 
grouping A are consistent with this grouping being a sub- 
cluster. 



4.25 EDCC 442 

From Figures lA25f a~l- (e) . EDCC 442 is an elongated cluster 
possessing a relatively dense core. Outside the central re- 
gions, there are five groupings with number densities greater 
than 40 gal/(fe~ 1 Mpc) 2 (labelled in Figure lAlST d^. Com- 
parison of Figures fA25f b~l and (c) shows that the cluster cen- 
tre of luminosity is closely correlated with the mean arith- 
metic centre due to the presence of several of the brightest 
galaxies apparently residing near or at the bottom of the 
cluster potential well. 

The ellipticity computation yields higher ellipticities at 
the larger mean radii due to the influence of the subgroup- 
ings, and smaller ellipticities in the central region, e.g., e = 
0.34 at R = 0.5 /i" 1 Mpc while e = 0.58 at R = 1.0 h' 1 Mpc. 
Note, however, that the position angle of the core galaxies is 
aligned with the subgrouping-influenced global cluster posi- 
tion angle of PA as 35°. Using the average ellipticities and 
position angles for the 0.5 and 0.75 /i -1 Mpc dispersion el- 
lipses to compute the core parameters of R c = 0.33 ft -1 Mpc 
and (Jo = 113±4 gal/(/i _1 Mpc) 2 yields a significantly better 
fit to the cluster density profile. 

The velocity distribution as shown in Figure lA"25lT ) is 
symmetric (S — < 0.01) but somewhat flattened compared 
to a normal distribution (K = 1.7 = 2. lax) with a slight 
indication of bimodality. Nevertheless, the Xvei test does not 
reject this distribution as normal (CL to reject the normal 
hypothesis = 51%). The three slowest galaxies with v < 
12, 500 km/s are probably in the foreground. 

Consistent with the appearance of the 2D visualization 
plots, the P test yields a statistically significant asymme- 
try of P/3 = 3% probability of random occurrence. The 
two-dimensional substructure reflects apparently real three- 
dimensional substructure in at least four of the subgroupings 
noted in Figure IA25f d'l. The nearest neighbor local velocity 
plot in Figure 1181 clearly shows that the galaxies in sub- 
groupings B, C, D, and E have local velocities consistent 
with their being members of bound subclusters. 

The kinematically-weighted centroid shifts tested by 
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the a algorithm do not indicate significant substructure 
(P a = 24.9%). However, with P K = 3.5%, the k test does in- 
dicate the probable presence of dynamically significant sub- 
structure. Correlating the bubble plot of Figure 

EUtf) with 

the segmentation/nearest neighbor plots of Figures lA25f d'l- 
(e) indicates again that the galaxies in groupings C and D 
are probable subclusters, and with the transition between D 
and E being gradual rather than sharp. Consistent with the 
"stringy" appearance of the local velocities of the nearest 
neighbors in Figure HS1 the bubble plot also gives some in- 
dication of core kinematics that may indicate non-virialized 
core dynamics. The lack of relatively large bubbles at the 
locations of groupings A, B, and E requires further inves- 
tigation to distinguish whether these are bound subclusters 
with velocity distributions indistinguishable from the global 
distribution, the local grouping motion is transverse to the 
line-of-sight, or the groupings are only a result of projec- 
tion effects (as appears likely for grouping A although this 
grouping does appear to contain a likely bound binary). 



5 SUMMARY OF RESULTS AND DISCUSSION 

To begin evaluating the physical significance of the results 
obtained here, consider that, in principle, any single test 
can detect the type of substructure to which it is sensitive 
regardless of the results of any other test. We believe that 
the methods used here to detect substructure in the galaxy 
distribution are generally accurate to the confidence levels 
stated due to the discriminating power of the 3-dimensional 
tests in substantially reducing (or eliminating) apparent 
substructure due to projection effects. In particular, the k 
test results combined with the bubble skyplot and plot of the 
local velocities/relative positions of the nearest neighbors 
appears to be a particularly useful combination for identify- 
ing true substructure. 

However, it should be remembered that it has been 
demonstrated through Monte Carlo simulation that sub- 
structure tests have both "failure to dete ct" and "false posi- 
tive" rates (see, e.g.. |Pinknev et alll996t) . Thus, before sum- 
marizing the results of the statistical tests applied to the 
present sample, it is appropriate to review estimates for the 
numerical reliability of the a and f3 tests. 

Pinkney et al. found that the a test is generally a ro- 
bust diagnostic of substructure, sensitive to dispersed sub- 
clusters, and that its effectiveness is not always diminished 
by a superposition of the core and a subcluster. In that 
study, the test identified substructure in 24 out of 36 cases 
of simulated mergers at a 10% probability of random oc- 
currence, second only to the A test. At a 5% probability 
level the test detected 16 out of 36 subclusters, and at 1%, 
12 out of 36. On the other hand, the test is also found to 
have relatively high false positive rates at the 10% and 5% 
significance levels. For example, in a cluster with 100 mem- 
bers, the test incorrectly detected substructure in 23 out 
of 100 cases at a 10% confidence level. At a 1% confidence 
level, however, its false positive rate w as closer to the other 
estimators tested (IPinknev et al.1 11996 - ). 

For the /3 test, the Pinkney et al. study found it to be 
the second best two-dimensional test out of the four two- 
dimensional tests that were evaluated; their most robust 2- 
dimensional test was the Fourier elongation test that was es- 



Table 6. Summary of statistical tests for substructure. The num- 
ber in the table is the percent probability that the test statistic 
is the result of a random occurrence. 
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sentially replaced in this study by the ellipticity algorithm. 
In that work, The f3 test identified substructure in 16 out 
of 36 cases with a 10% probability of random occurence, 
15 cases with a 5% probability, and nine cases with a 1% 
chance. The false positive rate was found to be 12%, 9%, and 
4% for significance levels of 10%, 5%, and 1% using a sample 
cluster size of 100 galaxies. Thus, it appears the test has a 
low detection rate and a high false positive rate. This was 
true for most of the one- dimensional and two-dimensional 
tests evaluated. As one would expect, the (3 test was effec- 
tive for detecting deviations from mirror symmetry, but not 
necessarily circular symmetry. In terms of merger scenarios, 
the test works best at early stages of merger viewed at right 
angles but not as well at later stages viewed nearly straight 
on. Perhaps most importantly, and not explicitly discussed 
in previous work, is the sensitivity of the results to the choice 
of cluster centroid. 

Given the above comments, it is obviously not possi- 
ble to make an absolute statement regarding the confidence 
level at which a particular test yields a true result. In many 
cases, a final judgment can be rendered by incorporating 
high-resolution, long-exposure X-ray data. In other cases, 
there may be no absolute determination possible. However, 
we believe that the approach used here to combine visualiza- 
tion plots with statistical tests and then to identify sources 
of positive signatures is robust in mitigating the relatively 
small probabilities for the statistical tests to yield a false 
positive or negative recognition. 

The results for the a, /3, k, and xtei tests are collected 
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Table 7. The number of clusters showing a substructure signal 
at less than or equal to a specified probability that the signal is 
a random occurrence. The entries include clusters whose Monte 
Carlo probabilities are rounded to the nearest percentile. 



Test 
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in Table |S| where it can be seen that only 4 of the 25 clusters 
failed to yield at least one signal at or below a 1% proba- 
bility for random occurrence (or, equivalently, at or greater 
than a 99% CL where the results have been rounded to the 
nearest percentile) . Table |7] tabulates the results in terms 
of the number of clusters that yielded substructure at less 
than the 10%, 5%, and 1% probabilities for random occur- 
rence. Multimodal central regions were detected in about 
half of the sample, and many of these appear consistent with 
merger scenarios. In general, the luminosity density in these 
clusters correlates closely with the number density, the few 
exceptions probably being due to the small level of incom- 
pleteness leaving out a few bright galaxies. However, this 
should be more closely examined in a future study. 

It is well known that giant elliptical galaxies in clusters 
usually reside in a region of high galaxy density, either in the 
central region or in bound subclusters, and it is believed that 
galaxy- environment interactions significantly affect galaxy 
formation and evolution. The absence of clear morphologi- 
cal data for most of the galaxies in this catalogue restricts 
analysis here to correlating brightness with dumpiness (as 
provided by the luminosity density contour plots). We find 
that 16 of the 25 clusters contain probable bound subclusters 
with luminosity density comparable to that in their cores. Of 
these, 14 clusters contain one or more subclusters in which 
the first, second, or third brightest galaxy is found. Espe- 
cially noteworthy are A957, A1663, and A1750 as clusters 
having subclusters comparable in size to the cluster core and 
containing a single extremely bright galaxy. The other 11 
clusters are A1238, A1620, A2814, A3027/APM268, A3094, 
S258, S301, S1043, APM933, and EDCC442. 

As mentioned in the Introduction, the clusters selected 
for this study from the 2dFGRS catalogue were required 
to have at least 70 members in order to have confidence 
in the statistical results, and to provide a small enough set 
to allow analysis of individual clusters. With the catalogue 
utilized for this study not containing clusters with member- 
ships exceeding N = 154, this selection criterion then nat- 
urally resulted in analysing low richness clusters of galaxies 
which are structures not as well studied as rich clusters and 
groups. Assessing the robustness of the approach here, it 
appears that smaller-sized groups in the 2dFGRS catalogue 
could reliably be analysed using the same methods. Future 
work with this particular set of clusters includes incorpo- 
rating d ata from the rece ntly published 2PIGG catalogue of 
clusters jEke et al.l l2004') and comparing against recent N- 
body simulations of cosmological concordance models. More 
detailed examinations are also planned of selected clusters 



by including new X-ray data, applying estimators to char- 
acterise cluster dynamics, and using analytical tools such as 
wavelet analysis or the KMM algorithm to further resolve 
smaller-scale structure. 

To summarize, we have presented a survey-level sub- 
structure analysis of 25 low richness clusters of galaxies con- 
tained in the 2dFGRS cluster catalogue. Most of the clusters 
possess features such as multi-component core/central re- 
gions, bound subclusters well-removed from the centre, pre- 
and post-merger indicators, and possible signatures of infall, 
rotation, or shear dynamics. Thus, we find that substructure 
in these clusters is the rule rather than the exception indi- 
cating that, as conjectured by others and supported by some 
N-body simulations, low richness clusters relax to structure- 
less equilibrium states on very long dynamical time scales 
(if at all). The results here also show that while doubling 
the number of redshifts available per cluster does not dras- 
tically alter the velocity statistics found in previous studies, 
the additional information affords a more comprehensive re- 
construction of cluster substructure. 
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velocity statistics and the k test, visualization of potential 
substructure correlated with the cluster velocity distribu- 
tion is accomplished through two additional types of plots: 
the peculiar velocity of each galaxy, v ga i — v c i , plotted in a 
3-dimensional position-velocity format, and the "local" ve- 
locity constructed from the average velocity of each galaxy 
and its nearest neighbors, vi oca i (= Vnn) — v c i, plotted in 
the same 3-dimensional position-velocity format. The zero 
point for the peculiar velocities is the cluster average veloc- 
ity, whereas the zero point for the local velocity is based on 
the local velocity of the galaxies residing within one core 
radius of the cluster centre. For these latter two plots, the 
position of each galaxy is marked by a filled circle when 
the k statistic for that galaxy has a < 5% probability of 
random occurrence and an open circle otherwise (the 5% 
probability threshold is arbitrary and can be varied to indi- 
cate the "strength" of a particular grouping). The vertical 
lines connecting the position of each filled circle galaxy to 
the midplane are for visualization purposes only. As seen in 
the individual cluster analyses, the local velocity plot pro- 
vides a useful visualization complement to the k statistic 
and bubble plot. Finally, plot (i) shows the histogram of the 
velocity distribution using a bin size of 150 km s _1 . 



APPENDIX A: VISUALIZATION PLOTS FOR 
THE CLUSTERS 

This appendix presents the set of visualization plots de- 
scribed in Section 3. For each cluster analyzed, plot (a) 
shows each galaxy position in (RA, Dec) converted to carte- 
sian (x, y) coordinates with superimposed dispersion ellipses 
for mean radii of 0.5 and 1.0 h^ 1 Mpc. The projected (2- 
dimensional) galaxy separations are inferred from the proper 
distance at photon time of emission, d p (t e ). Following the 
usual convention, north is along the positive y-axis and 
east is along the negative x-axis. The filled circles represent 
galaxies with velocities within 1.3<r of the mean velocity, 
whereas galaxies slower than 1.3cr are represented by a open 
circle, and galaxies faster than 1.3cr are shown with a plus 
sign. For a normal velocity distribution, approximately 10% 
of the galaxies should fall in the faster or slower group. Plots 
(b)-(e) are the 2-dimensional visualization plots described in 
Section 3.2. 

Plot (f) is the k "bubble skyplot" described in Section 
3.5 where the size of the circles centered on galaxy positions 
is proportional to the k statistic for each galaxy. As a com- 
plement to the substructure implications provided by the 
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Figure Al. Visualization plots for Abcll 930. 
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Figure A2. Visualization plots for Abcll 957. 
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Figure A3. Visualization plots for Abell 1139. 
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Figure A4. Visualization plots for Abell 1238. 



34 W. Burgett et al. 

(a) Galaxy Position 



1 



o 





-1 







— 




• 


+ + f. 

* 'Y^ 


\ **• 




• 


.-.."> 

••••_!» 


• 











-1 1 
Mpc 



(b) Number Density 



-1 




-1 



1 
Mpc 



(c) Luminosity Density 




o 





-1 



(d) Segmentation 











E 
• 

D 


A 
G 









-1 1 
Mpc 



(e) 4NN<0.25Mpc 



1 


-1 











E. .». 
* ^ '* 

- «uV * 


A 




t 


••• • (_ 













-1 1 
Mpc 



(f) k Bubble Skyplot 



-1 



o 



"ZD* 
D 



-1 



1 
Mpc 



(g) 3D Peculiar Velocity 



E 



CO 

O 




-2. 



(h) 3D Local Velocity 



(i) Velocity Distribution 




Velocity (1(f km/s) 



Figure A5. Visualization plots for Abell 1620. 
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Figure A6. Visualization plots for Abell 1663. 
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Figure A7. Visualization plots for Abell 1750. 
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Figure A8. Visualization plots for Abell 2734. 
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Figure A9. Visualization plots for Abell 2814. 
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Figure A10. Visualization plots for Abell 3027. 
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Figure All. Visualization plots for APM 268. 
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Figure A12. Visualization plots for Abell 3094. 
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Figure A13. Visualization plots for Abell 3880. 
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Figure A14. Visualization plots for Abcll 4012. 
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Figure A15. Visualization plots for Abell 4013. 



Substructure in 2dFGRS Clusters 45 



(a) Galaxy Position (b) Number Density (c) Luminosity Density 




-10 1 -10 1 -10 1 

Mpc Mpc Mpc 



(d) Segmentation (e) 4 NN < 0.25 Mpc (f) k Bubble Skyplot 




Figure A16. Visualization plots for Abcll 4038. 
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Figure A17. Visualization plots for Abell S141. 
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Figure A18. Visualization plots for Abell S258. 
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Figure A19. Visualization plots for Abell S301. 
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Figure A20. Visualization plots for Abell S333. 
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Figure A21. Visualization plots for Abell S1043. 
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Figure A22. Visualization plots for APM 917. 
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Figure A25. Visualization plots for EDCC 442. 



